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ABSTRACT 

An  analytical  and  numerical  study  is  made  of  a  low 
frequency  Sloan-Lawrence  structure  with  electrostatic  quadru*- 
poles  within  the  drift  tubes  for  accelerating  charged  micro¬ 
particles.  An  improved  set  of  equations  has  been  derived  to 
study  the  dynamics  of  the  charged  particles.  These  equations 
have  been  used  to  arrive  at  several  structure  designs,  and 
these  structures  in  turn  have  been  examined  for  angular  accep¬ 
tance,  phase  acceptance,  and  acceptance  in  charge  to  mass 
ratios.  The  ranges  of  operating  frequencies  and  charge  to 
mass  ratios  considered  are  30-50  KHz  and  26-84  coulombs/kg 
respectively . 
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CHAPTER  1 


INTRODUCTION 

This  thesis  deals  with  a  study  of  a  low  frequency 
Sloan-Lawrence  structure  for  the  acceleration  of  electrically 
charged  microparticles  for  the  simulation  of  micrometeoroids. 

The  structure  consists  of  a  series  of  cylindrical  drift  tubes 
and  accelerating  gaps.  In  this  type  of  structure  axial  par¬ 
ticle  stability  is  incompatible  with  simultaneous  radial 
particle  stability.  To  overcome  this  difficulty  and  to  provide 
stability  in  both  axial  and  radial  directions  the  effect  of 
auxiliary  focusing  elements  in  the  structure  is  examined. 

In  principle  there  are  various  kinds  of  focusing  arrange¬ 
ments  which  can  be  used,  such  as  solenoid  focusing,  quadrupole 
focusing  and  self  focusing.  Of  all  the  recently  developed 
focusing  systems,  quadrupole  focusing  is  the  most  efficient 
and  it  is  the  type  of  focusing  studied  in  this  work. 

The  idea  of  using  quadrupole  magnets  as  focusing  elements 
was  first  introduced  by  E.  D.  Courant  et  al^1^,  during  an 
investigation  of  the  feasibility  of  alternating  gradient 
synchrotrons.  In  19^1  D.  W.  Kerst  and  R.  Serber^^,  found  that 
in  the  induction  accelerator  the  strength  of  the  focusing 
force  is  limited  by  the  stability  criterion:  0<n<l,  where 
n  is  called  the  field  index  and  is  given  by  the  expression 
n  =  -^j  ,  r  is  the  radial  displacement  and  B  is  the  magnetic 

field.  If  n  is  much  greater  that  1  (rapid  decrease  of  field 
strength  in  the  radial  direction),  a  strong  focusing  field 


>  . 

' 

. 
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results  in  the  vertical  plane  (refer  to  fig.  1-1)  while  the 
radial  motion  is  strongly  defocused.  If  the  next  section 
has  a  large  negative  n  (rapid  increase  of  field  in  the  radial 
direction)  there  will  be  a  strong  radial  focusing  force  with  a 
large  defocusing  force  in  the  vertical  plane.  In  1951 
E.  D.  Courant  et  al  showed  that  the  net  effect  of  these  two 


Fig.  1-1  Section  of  Induction  Accelerator 

sections  is  convergent  in  both  the  vertical  and  the  horizontal 

planes  to  a  greater  extent  than  in  a  constant  gradient  magnetic 

field  with  a  value  of  n  lying  between  0  and  1. 

(9) 

In  the  same  year  J.  P.  Blewett  showed  that • the  same 
type  of  alternating  gradient  focusing  quadrupoles can  effec¬ 
tively  be  used  for  linear  particle  accelerators  of  the 
Alvarez  type.  His  lens  system  consisted  of  quadrupole  magnets 
placed  axially  inside  the  drift  tubes  of  the  accelerator  such 
that  in  alternate  drift  tubes  the  quadrupoles  are  rotated 
through  90  degrees  about  the  axis  of  the  accelerator.  The 
gradual  development  of  this  work  in  chronological  order  is 
listed  in  the  references  by  1,  3  to  12,  and  14. 
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This  work  consists  of  three  main  chapters.  The  first  of 

these.  Chapter  2,  deals  in  detail  with  the  equations  of  motion 

of  a  particle  in  the  accelerating  structure.  A  series  of 

equations  to  describe  the  particle  dynamics  is  derived.  These 

equations  are  more  accurate  and  general  than  those  used  by 

p  6  7  8  1  li 

previous  workers^5  ,J  J 

Chapter  3  consists  of  an  analytic  study  of  transverse 
particle  motion  in  a  long  accelerator  structure  under  the 
simplifying  assumptions  that  the  particle  velocity  remains 
constant  and  that  the  structure  is  periodic.  Analytical 
expressions  have  been  derived  for  the  transverse  excursions 
of  a  particle  at  the  injection  plane  and  at  all  succeeding 

i 

transverse  planes  at  the  middle  of  the  drift  tubes.  A  general 
expression  for  the  acceptance  in  phase  space  has  been  derived 
in  any  of  the  above  planes  along  the  accelerator. 

In  Chapter  4  a  numerical  study  of  the  behavior  of  various 
structures  is  made  using  the  equations  derived  in  Chapter  2. 
The  results  of  this  study  are  presented  in  tabular  and 
graphical  forms . 
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CHAPTER  2 


PARTICLE  DYNAMICS  IN  A  MICROPARTICLE  LINEAR  ACCELERATOR 

2.1  Radial  Motion  in  the  Accelerating  Gap 

2.1.1  General  Equations  and  Assumptions 

The  equations  of  motion  of  the  charged  microparticle  in 
the  accelerating  gap  are  derived  from  the  Lorentz  force 
equation : 


(2-1) 


where  t  is  the  time 

v  is  the  velocity  of  the  particle 
q  is  the  charge  on  the  particle 
B  is  the  magnetic  field  intensity 
E  is  the  electric  field  intensity. 

The  velocity  of  the  microparticles  is  very  low  and,  there¬ 
fore,  the  component  of  force  due  to  magnetic  field  is  neglected. 
It  then  follows  from  the  foregoing  equation  that  the  radial 
motion  of  the  particle  is  described  by 


mr-mr02  =  qEr(r,Z)  COS(cot  +  4>) 


(2-2) 


where  r  is  the  radial  excursion  of  the  particle  from  the  axis 
of  the  accelerator 

E^  is  the  radial  component  of  the  electric  field  intensity 


6  is  the  angular  velocity 
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z  is  the  distance  measured  along  the  axis  of  the  accel¬ 
erator 

m  is  the  operating  frequency 

$  is  the  phase  difference  between  the  instant  in  time 
at  which  the  accelerating  field  reaches  its  peak  value 
and  the  instant  at  which  the  particle  crosses  the 
electrical  centre(Z  =  0).  <J>  is  positive  if  the  particle 

crosses  the  electrical  centre  after  the  field  has 
reached  its  peak  value.  The  exact  definition  of  the 
electrical  centre  is  given  in  Appendix  1.  In  this 
work,  as  in  most  other  cases,  the  electrical  centre  is 
the  mid-point  of  the  gap. 

The  present  study  is  restricted  to  an  accelerator  of  the 
Sloan-Lawrence  type.  In  this  accelerator  the  reference  particle 
or  the  so-called  synchronous  particle  traverses  the  distance 
between  two  successive  gap  centres  in  exactly  T/2,  i.e.,  one 
half  cycle  of  the  applied  electric  field.  As  a  consequence 
the  synchronous  particle  traverses  all  the  gaps  while  the 
fields  are  accelerating  and  is  shielded  in  the  drift  tubes 
while  the  gap  fields  are  decelerating.  All  quantities  such 
as  energy,  velocity,  phase,  etc.,  associated  with  the  synchro¬ 
nous  particle  are  designated  by  the  subscript's'.  Fig.  2.1 
illustrates  the  principle  of  operation  of  a  Sloan-Lawrence 


structure . 
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phase 
=  0)t 


Fig.  2-1  Sloan-Lawrence  Structure 


In  computing  the  motion  of  a  particle  in  any  given 
accelerating  section  i.e.,  the  region  between  the  centres 
of  two  successive  drift  tubes,  the  following  simplifying 
assumptions  are  made: 

1.  It  is  assumed  that  the  particle  traverses  the  section 
at  constant  velocity.  This  is  justified  if  the  energy 
gained  by  the  particle  while  traversing  the  section  is 
small  compared  to  the  energy  which  the  particle  has 
when  entering  the  section. 


. 
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It  is  assumed  that  the  particle  traverses  the  section 
at  constant  radius.  This  is  justified  if  the  angle 
between  the  particle  trajectory  and  the  accelerator 
axis  is  small. 


>•  z 


z  = 


Note:  gn  is  the  length  of  n ^  gap;  LDn  is  the  length  of 

the  n^h  drift  tube;  LSn  is  the  length  of  n^  sec¬ 
tion;  A  and  B  are  the  mid  points  of  the  nth  and 
(n+l)tn  drift  tubes  respectively;  vs  n  is  the 
constant  velocity  of  the  synchronous ’particle  in 
a  single  section,  i.e.,  between  points  A  and  B. 

Fig.  2-2.  A  Single  Section  of  the  Sloan-Lawrence  Structure 


Fig.  2-2  shows  the  detailed  dimensions  of  a  single 
section  of  the  Sloan-Lawrence  structure.  The  synchronous 
particle  takes  a  half-cycle  to  go  from  one  gap  centre  to 
the  next.  As  a  consequence  of  this  and  the  assumption 
No.  1  stated  before,  the  particle  will  take  a  half-cycle 


8 


to  traverse  a  single  section.  Thus  the  length  of  the 

th  .... 
n  section  is  given  by 


LS  =  fv 

n  \  s  ,n/ 2 


(2-3) 


In  the  following  analysis  those  particles  which  are  not 
synchronous,  i.e.,  the  nonsynchronous  particles,  are 
treated  by  writing  their  velocity  vn  as 


v  =  v  /(1+k  )  - (2-4) 

n  s  ,n  n 

where  is  a  parameter  whose  value,  for  any  particle, 
changes  from  section  to  section. 

Using  2-4,  the  independent  variable  oat  of  2-2  can  be 
written  as 


<ot  -  (2n)z(i±!£B> 
V  T/  vs  ,n 


which,  by  use  of  equation  2-3,  becomes 


wt  =l§n)U+kn)  - (2-5) 


Because  of  the  assumption  2  of  the  previous  page  the  term 
containing  0  of  2-2  can  be  neglected. 

Thus,  substitution  of  the  relation  2-5  in  2-2  yields 


mr 


nz 


=  qEio(  r ,  Z  )  COS{^|  (l+kn)  +  <f>nJ. 


n 


n  '  n 
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Using  the  relation  r  =  vn  —  in  the  foregoing  equation 
one  obtains 


dr  = 


qEr ( r , Z ) 


mv 


n 


C0^Efn(1+kn)+‘t’n 


>dZ 


Integration  of  the  above  equation  from  -LSn/2  to  LSn/2 
(refer  to  Fig.  2.2)  yields 


LSn/2 


LSn/2 


Ar  = 

m  n 


dr 


= 


mv 


Er(r,Z)  COS- 


n 


nz 

LS 


( l  +  kn)  +  (J)n/  dZ 


ni 


n 


-LSn/2 


-LSn/ 2 


(2-6) 


where  Arn  is  the  difference  in  radial  velocity  of  the 
particle  at  the  beginning  and  the  end  of  the  n^h  section. 
From  A 1 —6  of  Appendix  1  it  is  seen  that  Er(r,Z)  is  an  odd 
function  of  Z.  Thus  2-6  becomes 


LSn/2 


Afn  =  - 


q  SIN  ( 4>n) 


mv 


Er(r,Z)  S!Nj (l+kn) |dZ 


nz 


n 


-LSn/2 


Assuming  kn<<l,  the  foregoing  equation  becomes 
q  SIN  ( 4>n) 


Arn  = 


mv 


n 


LSn/2 


Er(r,Z)  SIN(nZ/LSn)  dZ 


L-LSn/2 


+  k 


n 


LSn/2 

nz 

LSn 

-LSn/2 


Er(r,Z)  COS(nZ/LSn)  dZ 


(2-7) 


. 
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From  A2-3  and  A2-5  of  Appendix  2,  the  above  equation 
becomes 


Afn  =  - 


Vmq  SIN(<Dn) 


n' 


mv 


n 


R-knLSn 


dR 

dLSn 


(2-8) 


where  R  =  TSjn(0)  T1(nr/LSn) 

T  (r)  = 

s,n  ngn  i0(na/Lsn) 

2LSn 

Vm  =  peak  voltage  at  the  gap. 


Ig(IIr/LSn)  and  I1(IIr/LSn)  are  the  modified  Bessel 
functions  of  first  kind  of  zero  and  first  order  respec¬ 
tively  . 


the 


Tf  T  is 

nonsyn 

non synchronous 


defined  as  the  transit  time  factor 
particle  then  equation  2-8  becomes 


of 


Ar 


n 


Vmq  SIN(<f>n) 


mv 


n 


T 

nonsyn 


where 


T 

nonsyn 


R— kjqLS  ^ 


dR 

dLSn 


(2-9) 
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At  this  point  one  may  note  that  for  synchronous 
particle  or  almost  synchronous  particle  (kn-0),  the 
equation  2-8  becomes 


Ar 


n 


Vmq  SIN ( 8n ) 


mv 


n 


TSjn(°)  1 1  ( IIr/LSn) 


(2-10) 


Hence  it  is  seen  that  for  4>s  negative  Arn  is  positive, 
i.e.,  the  gap  field  is  defocusing.  On  the  otherhand,  if 
( ps  is  positive  Arn  is  negative  and  the  gap  field  is  focusing. 


2.1.2  Matrix  Representation  of  Particle  Motion 

The  entire  effect  of  the  accelerating  gap  on  the 
radial  motion  of  a  particle  is  assumed  to  take  place  at 
the  centre  of  the  gap  in  the  form  of  a  8-function.  If 
the  suffixes  In  and  2n  denote  quantities  just  to  the  left 
and  to  the  right  of  the  nth  gap  centre,  respectively,  then 
the  radial  motion  is  described  by 


r2n  =  rln 


(2-11) 


and 


Arn  =  r2n 


-  r 


In 


-Anr 


(2-12) 
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where,  from  2-8 


an  ■ 


Vmq  SIN(<t>n) 


mv 


n 


R-knLSn  dLgn 


(2-13) 


It  is  now  recalled,  that  for  the  purpose  of  computing 
the  radial  impulse  of  a  particle  it  was  assumed  that  the 
particle  traversed  the  section  at  constant  radius.  This 
constant  radius  is  taken  to  be  that  at  the  middle  of  the 
gap.  Hence,  it  follows  from  2-12 


Arn  =  -Anrp 


n 


(2-14) 


From  2-11  and  2-l4  the  radial  motion  of  the  micro¬ 
particle  at  the  gap  centre  can  be  written  in  the  form  of 
a  transfer  matrix  as 


(2-15 


2.2  Axial  Motion  in  the  Accelerating  Gap 
2.2.1  General  Equations  and  Assumptions 

The  assumptions  in  art. 2.1.1  will  also  be  used  in 
this  section.  The  equation  of  axial  motion  of  the  particle 
is  derived  from  the  Lorentz  force  equation  2-1  as 
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mZ  =  qEz(r,Z)  COS(u>t  +  <$> )  - (2-16) 

Prom  2-5  and  the  foregoing  equation  one  obtains  the 
equation  of  axial  motion  of  the  particle  in  the  n^ 
section  as 


mZ  =  qEz(r,Z)  COS 


nz 

LS 


( l  +  kn)  +  t})n 


(2-17) 


n 


Use  of  the  relation  Z  =  — 

dZ 


vn  in  the  above  equation 


leads  to 


dZ  = 


mv 


n 


/nz  ■ 

\LSn 


Ez(r,Z)  COS<  1-^-  (l+kn)+<|> 


'n 


dZ 


Integration  of  the  foregoing  equation  from  -LSn/2  to 
LSn/2  yields 


LSn/2 


-LSn/2 


LSn/2 


q 


mv 


Ez(r,Z)  COS 


n 


nz 


LS 


( l+kn )  +(j)n  dZ 


n 


-LSn/2 


(2-18) 


where  AZn  is  the  difference  in  axial  velocity  of  the  par¬ 
ticle  at  the  beginning  and  the  end  of  the  n^  section. 
Using  Al-5  of  the  Appendix  1  the  foregoing  equation 


become  s 
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AZn  = 


LSn/2 


LSn/2 


Ez(r,Z)  COS^  dZ  -  k 


LS 


n 


Ln 


-LSn/2 


f  [§}E^>Z)  SIN(§J 


-LSn/2 


(2-19) 


Integration  (refer  to  A2-7  and  A2-8  of  Appendix  2) 
of  the  above  equation  yields 


Vmq  COS  ( <j)n ) 
mvn 


n 


(r)  +  knLS 


dTs ,n(r) 
n  dLSn 


(2-20) 


In  all  numerical  calculations  the  above  change  in 
velocity  is  assumed  to  take  place  at  the  middle  of 

the  gap  in  the  form  of  a  6-function. 


2.2.2  Energy  Gain  of  a  Microparticle 


The  energy  gain  AEn  of  a  microparticle  in  traversing 


the  n^*1  section  is 


AE„  =  q 


n 


LSn/2 


Ez(r,Z)  COsjj!|-  (l+kn)+cf>n  j  dZ 


-LSn/2 


+  /  Er(r,Z)  cos|~  ( l  +  kn)  +cf)n 
1 


dr 


dZ 


(2-21) 


V 
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where  r-^  and  r^  are  the  radial  excursions  of  the  particle 
at  the  beginning  and  the  end  of  the  n^  section  respectively. 
If  r  is  the  first  order  derivative  of  r  w.r.t.  Z  then 

dr  =  v  dZ 


Using  this  relation  in  2-21  and  assuming  v  is  constant 
in  the  section,  one  obtains 


AEn 


i  fnz 

pr 


(l  +  kn)  +  <J)n  \  dZ 


+  r  /  Er  ( r ,  Z  )  COSJ^  (l+kn)  +  <|> 


Jnz 


-LSn/2 


1LS 


n'  Yn/ 


n 


dZ 


Note  that  the  radial  force  is  computed  as  previously, 
i.e.,  by  assuming  that  the  particle  traverses  the  gap  at 
a  fixed  radius  r.  However,  to  compute  the  energy  gain, 
the  particle  is  permitted  to  undergo  a  radial  displace¬ 
ment  of  r/  dZ .  Upon  integration  of  the  above  two  integrals 
(refer  to  2-18  and  2-20  for  the  first  integral  and  2-6, 

2-8  for  the  second  integral)  the  above  equation  becomes 


■ 
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AEn  =  vm<3 


COS  ( 4>n) 


dT„  n(r) 

Ts  n(r)+knLSn— - 

n  dLSn 


-  r'  SIN(*nK  H-knLSn  dR 


dLS 


n 


- - - (2-22) 


2.2.3  Incompatibility  between  Radial  and  Axial  Stability 

The  condition  of  radial  stability,  i.e.,  the  radial 
focusing  of  an  accelerating  particle  has  already  been 
discussed  in  section  2.1.1.  The  condition  of  axial 
stability,  i.e.,  the  presence  of  net  axial  restoring 
force  toward  the  position  of  the  synchronous  particle  is 
derived  as  follows: 

The  Z  coordinate  of  the  synchronous  particle  may  be 
written  as 


Z  =  Z  g  +  6 


where  6  is  the  Z  coordinate  of  the  nonsynchronous  particle 
referred  to  that  of  the  synchronous  particle. 

Equation  2-17  then  becomes 


m( Z  s+6  ^ 


qEz(r,Zs+6 ) 


COS- 


n(zs+6) 


LS 


n 


( l+kn )  +4> 


Expanding  Ez(r,Zs+6)  in  Taylor’s  series  about  Ez(r,Zs)  and 


neglecting  the  higher  order  terms  of  6  the  above  equation 
becomes 


•  •  •• 


m( Zg+6  ) 


=  q 


Ez  (r  ,ZS  )  +  <S 


3EZ  ( r ,  z  s  ) 


BZ 


COS-^ 


n(z  +6) 


LS 


-+<f> 


n 


n 


,  x  ncz.+s) 

=  qEz(r,Zs)  COsJ__s - +4, 


LS 


n 


n 


3E7(r,Z  ) 

+  q6 — - cos^ 

az  o 


n(z  +6) 

s +4> 


LS 


n 


n 


If  one  assumes  AZ<<ZS3  kn<<l,  and  4>n-c})Sj  then  the 
foregoing  equation  can  be  written  as 


3E  7  ( r ,  Z  o  ) 

m6  =  q6 - —  COS 

azQ 


nz, 


LS 


-+4>: 


n 


Use  of  the  relation  6  =  ^  v  v.  in  the  above 


equation  yields 


dZ  s  s»n 


q6  3E7(r,Z  ) 
d6  = - z - —  C0S4 


mv 


s ,  n 


az 


nz 


LS 


-+<$>sy  dz 


n 


Integrating  the  above  equation  from  -LSn/2  to  LSn/2, 
one  obtains 


LSn/2 


LSn/2 


A6n  = 


d6  = 


8EZ ( r ,  Zs ) 


mv 


s ,  n 


az 


nz„ 

COS^ — -+<f)s>  dz 


LS 


n 
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where  A6n  is  the  difference  in  relative  axial  velocity 
of  the  particle  with  respect  to  synchronous  particle  at 
the  beginning  and  the  end  of  the  n^h  section.  Assuming 
that  the  particle  crosses  the  section  at  constant  6  and 
using  the  result 


0 


the  foregoing  equation  becomes 


LSn/2 


q5  SIN((f)„) 


mv 


s  5n 


9EZ(r,Zs) 


-LSn/2 


s 


Upon  integration  by  parts  and  assuming  Ez  is  zero  at 


_LSn/2 ,  the  above  equation  can  be  written  as 


LSn/2 


A5 


mv 


s  ,n 


nz. 


SIN(jis)^n  j  j  Ez  ( r ,  Zg  )  COSii|s.  dZ 


LS 


n 


-LSn/2 


Use  of  A2-7  of  Appendix  2  in  the  above  equation 
yields 


A6 


n 


nq6  SIN (  4>s  ) 
LSn- mvs 


/ 

<  T 


(r )  V 
s,m  '  m 


(2-23 


» 


19 


From  the  foregoing  equation  it  is  seen  that  for  4>s 
positive  A6n  is  positive,  i.e.,  the  gap  field  forces  are 
defocusing  and  the  particle  under  consideration  moves  away 
from  the  synchronous  particle.  On  the  other  hand  for  <J>S 
negative  A6n  is  negative,  and  the  above  effect  is  reversed. 
From  the  foregoing  result  and  by  recalling  that  the  radial 
focusing  occurs  for  positive  4>s  it  is  concluded  that  radial 
stability  and  axial  stability  cannot  be  achieved  simulta¬ 
neously  for  the  type  of  accelerating  gap  considered.  This 
phenomenon  of  incompatibility  is  known  as  Macmillan’s 
theorem.  As  a  compromise  between  those  two  conflicting 
stability  conditions  the  synchronous  phase  <j>s  is  always 
chosen  close  to  zero.  There  are  no  suitable  methods  which 
can  eliminate  the  axial  instability  of  particles  if  cj)s  is 
positive.  However,  there  are  various  methods  which  can 
be  used  to  eliminate  the  radial  defocusing  effect  of  the 
gap  if  <J>S  is  negative.  For  this  reason  4>s  is  always' 
chosen  negative.  In  that  case  small  radial  defocusing 
of  the  gap  field  can  be  compensated  for  by  an  auxiliary 
focusing  system.  In  this  study  the  auxiliary  system 
consists  of  electrostatic  quadrupoles  which  are  placed 
co-axially  within  the  drift  tubes. 


-  I  £  I 
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2.3  Motion  of  a  Particle  in  a  Drift  Tube  with  an  Electro¬ 
static  Quadrupole 

An  electrostatic  quadrupole  consists  of  4  rectangular 
hyperbolic  electrodes  which  are  alternately  positive  and 
negative,  and  are  symmetrically  placed  about  the  axis  of 
the  accelerator  as  shown  in  Fig.  2-3. 


y 


>-  x 


Fig.  2-3  Sectional  View  of  the  Field  Lines  Between 


Quadrupoles 


A  solution  to  Laplace’s  equation  in  the  region  between 
the  hyperbolic  surfaces,  and  which  can  be  made  to  satisfy 
the  appropriate  boundary  conditions  is  given  by 


f  =  |k  (x2-y2) 


( 2-24 ) 


■ 
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Since  ¥  =  +V  at  —  -  =  1 

b2  b2 

2  2 

and  Y  =  -V  at  2L_  -  Z_  =  _>i 

b2  b2 

one  obtains  that 

k  =  iV - (2-25) 

b2 

where  2V  is  the  voltage  between  the  quadrupoles. 

Hence,  from  E  =  -VT  it  follows  that 

Ex  =  -  2VX  - (2-26) 

b^ 

Ey  =  ^|y - (2-27) 

b2 

Consequently,  the  equations  of  motion  for  a  particle  of 
charge  q  are 

x  =  -Q2x 

- (2-28) 

y  =  Q2y 

Q2  =  2V  a  - 

b2  m 


where 


(2-29) 


' 

- 
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Evidently,  for  a  positively  charged  particle  the 
motion  in  the  x-direction  is  focused,  while  the  motion  in 


the  y-direction  is  defocused. 

The  solutions  to  2-28  are 

x  =  a-L  COS  ( Qt )  +  a2  SIN(Qt)  - (2-30) 

y  =  a3  COSH(Qt)+a4  SINH(Qt)  - (2-31) 

and  x  =  -a1Q  SIN(Qt)+a2Q  COS(Qt)  - (2-32) 

y  =  a3Q  SINH(Qt)+a4Q  COSH(Qt) - (2-33) 


where  a2,  a2,  a3  and  a4  are  constants. 

With  the  assumption  that  the  velocity  vn  of  the 
particle  remains  constant  while  passing  the  (n^)  quadru- 
pole,  one  may  write  the  time  spent  in  the  (n^)  quadrupole 
as 


where  LQn  is  the  length  of  the  nth  quadrupole. 

If  t  =  0  when  the  particle  enters  the  quadrupole, 
then  at  t  =  Tn/2  the  particle  will  have  reached  the 
transverse  plane  at  the  middle  of  the  quadrupole. 


- 
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Denoting  the  quantities  at  the  entrance  and  the  midplane 
of  the  quadrupole  by  subscripts  0  and  1  respectively,  one 
obtains  from  2-30,  2-31,  2-32,  and  2-33 


Xq  a 2  ,  Xq  -  a2<3 


y0  a3’  y0  a4Q 


/QTn\ 

aq  COS  [~^J  +  a 2 


x 


1 


SIN 


'QT 


n 


X1  = 


-a^Q  SIN 


Pr)  +  a2Q  cos&) 


y1  =  a3  COSh(5!s)  +  a„  sinh(3Pi) 


yl  ~  a^Q  SINH\^  -  J  +  a^Q 


COSH 


© 


Prom  the  above  relations  one  finally  obtains  the 
motion  of  the  particle  in  the  first  half-length  of  the 
quadrupole,  in  matrix  form  as  follows: 
in  the  focusing  plane 


( 2-34 ) 


' 

. 
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in  the  defocusing  plane 


yx\  /  cosh 


QT 


n 


2 


l'n 


\  ^n 

\Q  SINH  - 


1  SINH  ^  n\ 
Q  2 

1*0) 

fyo\ 

QT 

■  td 

COSH  —J 

w 

n 

(2-35) 


where  T  and  are  called  the  transfer  matrices  in  the 
focusing  and  the  defocusing  planes  respectively. 

It  can  readily  be  shown  that  if  T  ^  and  are  the 

transfer  matrices  in  the  focusing  and  the  defocusing  planes, 
respectively,  of  the  second  half-length  of  the  quadrupoles 
then 


T  1  =  T  - - - (2-36) 

F  F 

T  /  =  Td  - (2-37) 


2.4  Transfer  Matrix  for  the  Radial  Motion  of  a  Particle 
for  One  Repeat  Section 

In  this  study  the  quadrupoles  are  arranged  so  that 
the  axes  x-y  are  rotated  through  90°  in  successive  drift 
tubes.  Thus  the  particle  traveling  in  the  focusing  plane 
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(y=0  plane)  in  one  drift  tube  will  travel  in  the  defocusing 
plane  (x=0  plane)  in  the  next  (refer  to  Pig.  2-4)  so  that 
the  focusing  and  the  defocusing  action  in  a  plane  is 


x 


h 


th  ,  . 
n  section 


y 

A 


\  th 


(n+1; 
section 


y 


7^ 


-4 


Z 


:z 


1  2 


5  6 


1  x  I 

3  4 

n^h  Repeat  Section 


I  y*T  I 

7  8  9 


"Quadrupoles 


Fig.  2-4  One  Repeat  Length  of  the  Quadrupole 
Focused  Linear  Accelerator 


repeated  after  alternate  drift  tubes.  The  length  between 

alternate  drift  tubes  is  known  as  one  repeat  length.  The 

transfer  matrices  T  ,  T  ,  TJ  ,  T  ,  T  ,  and  T  for  the 

regions  1-2,  4-5,  5-6,  8-9  and  the  gap  centres  3  and  7 

respectively  have  already  been  derived  and  are  shown  in 

equations  2-34,  2-35  and  2-15-  To  formulate  the  transfer 

matrix  of  the  complete  repeat  section  shown  in  Fig.  2-4 

23  34 

one  has  yet  to  find  out  the  transfer  matrices  T  ,  T  , 

T^  and  T^  of  the  drift  regions  2-3,  3-4,  6-7  and  7-8 
respectively.  Since  the  particle  traverses  these  drift 
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regions  at  constant  velocities  one  can  write  the  equation 
of  radial  motion  for  a  particle  in  the  region  2-3  as 


•  o 

r3  «  r2 

r3  =  r2  +  r2t23 


where  t23  =  L23/Vn  and  L23  is  the  length  between  the  sec¬ 
tions  2  and  3  in  Fig.  2-4. 

The  foregoing  equations  are  written  in  the  matrix 
form  as 


(2-38) 


Thus  the  transfer  matrix  T23  for  the  radial  motion  of  the 
particle  between  the  sections  2  and  3  is 


/i  t23\ 

\0  1  / 


(2-39) 


In  a  similar  way  one  can  find  the  transfer  matrices 
T^\  T^  and  T^ .  It  is  evident  that  the  above  transfer 
matrix  2-39  for  r  and  r  is  also  valid  for  x  and  x,  and 
y  and  y.  Hence,  the  equation  describing  the  motion  parallel 


' 
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to  the  plane  of  the  paper  In  Pig.  2-4,  i.e.,  in  the  verti¬ 
cal  plane  and  between  the  sections  1  and  5  can  be  written 
as 


/*i\ 

t45  t31)  t3  t23  t12 


- (2-1(0) 


where 


T 


aN 


t45  t34  t3  t23  t12 


(2-41) 


Similarly  the  matrix  equation  between  the  sections 
5  and  9  is 


=  t89  t78  t7  t67  t56 

1  y5  ' 

=  T 

ly5  \ 

BN 

w 

N 

\  K ! 

\  5  /j 

(2-42) 


where 

T  =  t89  t78  t7  t67  t56 
BN 


(2-43) 
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From  2-40  and  2-42  the  equation  of  motion  between  the 
sections  1  and  9  in  the  vertical  plane  is 


x9  ' 


'  x  9  N 


=  T  T 
gN  aN 


I 


'  1  'N 


/  x,  \ 


=  T. 


'N 


1 


\X1 


1  'N 


(2-44) 


where 


T  =  T  T 
N  $N  aN 


(2-45) 


T  is  the  transfer  matrix  for  the  n"^  repeat  section 
in  the  vertical  plane. 

It  is  evident  from  the  foregoing  analysis  that  the 
particle  motion  perpendicular  to  the  plane  of  the  paper 
in  Fig.  2-4,  i.e.,  in  the  horizontal  plane  and  between 
the  sections  1  and  9  is  given  by 


/  \ 


V  9/N 


/y 


=  T 


'N 


1 


(2-46) 


\yi/N 


where  T  =  T  „T  T 
N  aN  3N 


/  t  h 

T„'  is  the  transfer  matrix  for  the  N 


(2-47) 


in  the  horizontal  plane. 


repeat  section 
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CHAPTER  3 


SOME  PROPERTIES  OF  PARTICLE  MOTION  IN  THE  TRANSVERSE  PLANE 

IN  A  LINEAR  ACCELERATOR  WITH  ELECTROSTATIC  QUADRUPOLES 

3.1  General  Equations 

and  ,  the  transfer  matrices  for  the  transverse 
motion,  change  from  one  repeat  section  to  the  next.  So 
it  is  very  complicated  to  study  the  transverse  behaviour 
of  the  beam  analytically.  In  the  next  chapter  a  numerical 
study  of  the  beam  properties  will  be  made,  taking  into 
proper  account  the  variation  of  T^  and  T^  with  N.  In  the 
present  chapter  the  beam  is  studied  analytically  with  the 
simplifying  assumption  that  T^  and  T^  are  the  same  for 
all  N.  This  approach  does  not  yield  an  exact  answer, 
however  it  provides  a  great  amount  of  physical  insight 
into  the  transverse  behaviour  of  the  particle  beam.  In 
particular  then  the  assumptions  made  in  this  chapter  are 
the  following: 

1.  The  particle  traverses  the  accelerator  at  constant 
velocity  and  phase. 

2.  The  repeat  sections  are  identical  in  all  respects  and 
symmetric  about  the  middle  plane  transverse  to  the 
axis  of  the  accelerator. 

By  the  foregoing  assumptions  one  can  write  T^  and 

V  as 


' 


■ 
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(3-D 


/ 

T 


(3-2) 


where  T  ,  T  ,  T21  and  T22  are  the  elements  the 
transfer  matrix  T  of  one  repeat  section  between  focusing 
and  focusing  and  T^,  3  ^21  an<^  ^22  are  the  elements 

of  the  transfer  matrix  T/  of  one  repeat  section  between 
defocusing  and  defocusing. 

Now  if  the  vertical  planes  in  the  first  drift  tube  are 
assumed  to  be  focusing  then  the  vertical  planes  in  each 
alternate  drift  tube  will  also  be  focusing.  Thus  one  can 
write 


where 


/x 

lx\ 

and 

\xJ 

n  W 

are  the  coordinate  vectors  in  the 


P 
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vertical  planes  and  at  the  beginning  of  the  first  and  the 
t  h 

end  of  the  p  repeat  section  respectively. 

By  use  of  3-1  the  foregoing  equation  becomes 


/  x  \ 


(3-3) 


Similarly  in  the  horizontal  and  the  defocusing  planes 


jy\ 

y  \ 

=  T/p 

W 

P 

y  L 

(3-4) 


where 


and 


are  the 


coordinate  vectors 


in  the 


horizontal  planes  and  at  the  beginning  of  the  first  and 
the  end  of  the  p^*1  repeat  section  respectively. 


3.2  Radial  Stability 

By  definition  radial  stability  means  that  all  the 
elements  of  the  transfer  matrices  TP  and  T/  P  remain  bounded 
as  p  increases  indefinitely.  This  condition  is  necessary 


and  sufficient . 
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The  required  condition  for  stability  for  is  arrived 
at  in  the  following  way: 

If  A  is  one  of  the  eigenvalues  of  the  matrix  T,  and  if  use 

13 

is  made  of  the  relation  det  T  =  1,  then  the  characteristic 
equation  of  T  is 


A  -ACT  +T  ) 


+  1=0 


or , 


(T  +T  )+ 
K  11  122)~ 


/ 


(Tll+T22 


(3-5) 


where  i  =  1,2 

To  diagonalize  T  it  is  written  as  follows 


T  =  ADA 


(3-6) 


where  A  is  the  transform  matrix  and  D  is  the  diagonal 

IX  0 

matrix  with  Aj_  as  the  diagonal  elements  i .  e .  ,  D  = 

\0  A  2 


Hence 


DP 


(3-7) 


where  p  is  an  integer. 
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Thus  from  3-6  and  3-7 


Tp  =  A(AP)  A-1  - (3-8) 


P  r, 

Thus,  A ^  is  the  only  p  dependent  term  in  the  matrix  T  . 
This  means  that  the  elements  of  the  transfer  matrix  T^ 
will  be  finite  for  any  p  if  and  only  if 


A.|  <  l - (3-9) 


Thus  from  3-5  and  3-9  one  obtains 


Tll+T22 


±  / (Tl1^22)2-: 


<  1 


(3-10) 


It  follows  that 


T  +T 
ill 22 


<  2 


(3-11) 


implies  stability, 


t11+t22 


2 - (3-12) 


implies  instability. 


■ 


Similarly 


3^ 


(3-13) 


assures  stability  and 


/ 


(3-14) 


means  instability. 

3.3  Equation  of  Transverse  Motion  for  a  Stable  Particle 

Because  of  the  assumptions  of  periodicity  in  section 
3.1,  and  T  of  equations  2-4l  and  2-43  respectively 

can  be  written  as 


T 


al 


T 


T 


T 


a 


(3-15) 


and 


/$H  312\ 


T  =  T 


T  =  T 
3N  3 


-(3-16) 
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\32i  3  2  2 / 


35 


where  ai2>  a21  anci  a22  are  t^ie  elements  of  the  trans¬ 

fer  matrix  Ta  of  one  half  repeat  section  between  focusing 

to  defocusing.  3pp,  ^12*  ^21  an^  ^22  are  elemen"1:s  0 f* 

the  transfer  matrix  of  one  half  repeat  section  between 
defocusing  to  focusing. 

It  can  be  shown  that  since  the  repeat  section  is 
symmetric  about  the  middle  plane  transverse  to  the  accele- 


rator 

axis  the  following 

i — 1 
i — 1 

3 

“12\ 

22 

^cvP" 

i — i 

oa 

i — 1 

C\J 

a22/ 

i — i 
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From  3-18  and  3-19  it  follows  that 

T]_]_  =  T22  =  Tii  =  T22  =  alla22+a12a21 - (3-20) 

T12  =  a22a12  _ _ ,  ft) 

T21  alla21 
T12  alla12 

-7-  =  - - (3-22) 

T21  a22a21 

From  the  stability  conditions  3-11,  3-13  and  the  equation 
3-20  one  can  write  for  a  stable  particle 

Ti  -]  +T77  T11+T77 

— — —  =  — - —  =  COS(y)  - (3-23) 

2  2 

where  y  is  real. 


From  3-20  and  3-23  it  follows  that 


T11  =  T22  =  T11  =  T22  =  cos(t)  — — - (3-24) 

From  3-24  and  the  fact  that  det  T  =  det  T/  =  l^1^)  0ne 
can  write 

t12t21  =  "  SIN2  ( y  )  - - - -—(3-25) 
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T12T21  =  -SIN2(m) 


(3-26) 


To  find  T^2  and  T 2 ,  3-21  is  written  as 


hi  =  a22a12  =  _g2  - ( 3-27 ) 

T21  alla21 


From  3-25  and  3-27  one  obtains 


T12  =  3  SIN  ( y  )  - - - (3-28) 

T  =  SIN ( y  )  - - - (3-29) 

21  6 


From  3-22  and  3-27  it  follows  that 


(3-30) 


If  one  takes 


a22  =  02 
all 


(3-31) 


equation  3-30  becomes 
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(3-32) 
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Then  substitution  of  3-32  in  3-26  yields 

=  i-  SIN(u)  - 

9  ^ 

t21  *  -  J-  SIN(p)  - 


(3-33) 

( 3-3*0 


Thus  from  3-24,  3-28,  3-29,  3-33  and  3-34  the  transfer 
matrices  T  and  T/  can  be  written  as 


COS ( y )  6  SIN ( y  )\ 


T  = 


V--  SIN(y) 

e 


cos 


(y) ) 


(3-35) 


y 

T  = 


COS(u)  \  SIN(u) 


(3-36) 


■f-  SIN(y) 

P 


COS(y) 


Thus,  from  3-3  and  3-4  and  the  foregoing  equations  the 
equations  of  transverse  motion  of  the  stable  particle  are 
x  \  / x  v  /  COS(py)  BSIN(py)\/x 

=  Tp  I  1=1  II  I - (3-37) 
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n  >-i  siN(pu)  cos(py)/\x/0 

U  P 
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ly\ 
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—  rpP 
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\yl 
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/X 

/  y\ 

COS(py) 


fl2 

SIN (py  ) 


- (3-38) 


where 


and 


\x/o 


y 


are  the  coordinate  vectors  of  the 


0 


respectively  and 


injected  particle  in  the  vertical  and  the  horizontal  planes 

x\  /y  \ 

and  are  the  coordinate  vectors 

*4  'y/p 

of  the  particles  at  the  end  of  the  p^h  rep6at  section  and 
in  the  vertical  and  the  horizontal  planes  respectively. 

Similarly  the  equations  of  transverse  motion  at  the 
middle  of  any  repeat  section  are 
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- (3-39) 
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(3-40) 
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fy  \ 

!*\ 

where 

1  and 

. 

u 

middle  of  the  (p+l)th 
and  horizontal  planes 


are  the  coordinate  vectors  at  the 

repeat  section  and  in  the  vertical 
respectively . 


3.4  Acceptance  of  Particles  in  the  Transverse  Plane 

In  order  that  the  particles  are  transmitted  through 
the  length  of  the  accelerator,  it  is  necessary  that  the 
finite  excursions  of  the  particles,  as  described  by  equa¬ 
tions  3-37  and  3-38,  do  not  exceed  the  distance  "b"  from 
the  axis  of  the  accelerator  to  the  tips  of  the  quadrupole 
electrodes.  The  object  of  this  section  is  to  compute  the 
initial  conditions  of  the  particles  so  that  the  above 
boundary  condition  is  not  violated. 

By  elimination  of  py  from  equations  3-37  and  3-38  of 
the  previous  section  one  obtains 


Xp2  +  B2ip2  =  X02  +  62x02 


(3-41) 


2  8*2  2.8*2 
yP  +  yr  yp  =  yo  +  pr  yo 


(3-42) 


Similarly  by  eliminating  py  from  equations  3-39  and 


3— 4 0  one  obtains 
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yp+h  (a22+e2a21)  +  (a12  +  B\Pyp+k  ■2(a12a22+g2allcl21)yp+)sV!S 


2.,  2  \  •  2 


=  XQ  +e2  io  - (3-43) 

and 

2^2  2  2 

xp+!s  (all+~4c‘21)+xp+%  (a12+V22)_2xp+35xp+»5(alla12+V21a22) 
0  0  ^ 


2  .  B2-2 

yo  +  Vo 

04 


(3-44) 


The  equations  3-41  and  3-42  represent  the  equations 
of  ellipses  whose  major  and  minor  axes  are  coinciding  with 
the  coordinate  axes,  and  the  equations  3-43  and  3-44 
represent  the  general  equations  of  ellipses  whose  axes 
rotated  clockwise  through  the  angles  4>x  and  y  respectively 
where 


\  tan 
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By  use  of  the  relations  3-27  and  3-31,  and  the  relation 
det  T~  =  1  the  equations  3-43  and  3-44  become 
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(3-47 ) 
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(3-48) 
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and 


<t  =  90°,  *  =  90° 

x  y 

Now  from  3-41  and  3-42  it  Is  seen  that  corresponding 
to  each  set  of  (xp,  xp )  and  (yp>yp)  there  corresponds  a 
whole  family  of  sets  (xp,xp)  and(yQ}yp).  The  limiting 
two  families  of  sets  (xq,Xq)  and  (yg,;^)  correspond  to  the 
boundary  conditions  (xpmax,0)  and(ypmax,0)  and  are  conse¬ 
quently  given  by  the  relation 
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The  above  two  equations  represent  two  ellipses  in  the 
phase  space  (Xq,Xq)  and  (y0,y0)  respectively.  Any  particle 
whose  initial  (xq,Xq)  coordinates  lie  within  the  ellipse 
defined  by  the  equation  3-49  will  be  transmitted  through 
the  accelerator  such  that  the  x-coordinates  at  the  end  of 
any  repeat  section  will  always  satisfy 


x  <x 

pmax 

Similarly,  any  particle  whose  initial  (yQ,^) 
coordinates  lie  within  the  ellipse  defined  by  equation 
3-50  will  be  transmitted  through  the  accelerator  such 
that  the  particle's  y-coordinate  at  the  end  of  any  repeat 
section  will  always  satisfy 


y  <y 

J  Jpmax 


Similarly,  from  3-47  and  3-48  one  obtains 


2 

x0 


•  2 

x0 


(  0  y  ,  ^  ) 

^p+^max 


2 


e  y 


2 


=  1 


p  +  igmax 


$ 


-(3-51) 


and 
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2  •  2 

yp _ +  yp 


/xp+i$max\  2 

xp  +  ismax\ 

\  0  ' 

l  3  ) 

i 


(3-52) 


where  xp+^max  and  yp+^max  are  the  maximum  excursions  at  the 
middle  of  any  repeat  section  and  in  the  horizontal  and  the 
vertical  planes  respectively.  Here  again  the  two  ellipses 
in  the  (xq,Xq)  and  (yQ,^)  planes  define  the  limiting 
injection  parameters  such  that  no  particles  exceed  xp+pmax 
and  yp+^max  at  the  middle  of  any  repeat  section.  Now,  it 
is  customary  to  define  the  acceptance  of  an  accelerator 
in  any  given  phase  space,  for  instance  Xq,Xq,  as  1/IT  times 
the  area  of  the  corresponding  limiting  ellipse. 

Thus,  the  acceptance  at  the  entrance  to  the  accelera¬ 
tor  and  in  the  vertical  plane  is 
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Similarly,  the  acceptance  at  the  accelerator  entrance 


and  in  the  horizontal  plane  is 
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Prom  3-53  and  3-54  one  obtains 
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pmax 


0  yp+i^max 


(3-55) 


and 
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pmax 


xp+^max 

0 


(3-56) 


Now  from  Fig.  2-4  of  section  2-4  it  may  be  seen  that 
the  vertical  motion  of  a  particle  will  take  place  alter¬ 
nately  in  focusing  and  defocusing  planes.  Similarly  the 
horizontal  motion  will  take  place  alternately  in  the 
defocusing  and  the  focusing  planes.  It  is  evident  from 
Liouville's  Theorem  that  the  phase  space  acceptance  for 
vertical  motion  at  the  entrance  of  the  accelerator  and 
in  the  focusing  plane  is  equal  in  magnitude  to  the  phase 
space  acceptance  for  vertical  motion  one  half  repeat 
section  later  and  in  the  defocusing  plane.  A  similar 
conservation  of  phase  space  acceptance  holds  for  motion 
in  the  horizontal  plane.  In  view  of  this  and  because  the 


I 
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structure  is  assumed  to  be  periodic,  it  follows  that  the 
acceptance  at  the  entrance  to  the  accelerator  has  the 
same  magnitude  for  horizontal  and  vertical  motion  respec¬ 
tively  ,  i  .  e  .  , 


(3-57) 


Thus,  using  the  equations  3-53  and  3-54  in  3-57  one 
obtains 


xpmax  ®  yp+^max  ®  ^pmax  xp+^max 


(3-58) 


Now  using  the  boundary  condition  that 


xpmax  b 


and  the  relation  3-58  the  equations  3-49,  3-50,  3-51,  and 
3-52  become 
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The  ellipses  described  by  the  equations  3-59  and  3-60 
are  shown  in  Fig.  3-1  below 
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x0>y 


0 


Fig. 


3-1  Transverse  Acceptance  at  the  Plane  of  Injec¬ 
tion  of  an  Accelerator  with  Identical  Repeat  Sections 


Now,  a  cylindrically  symmetric  beam  is  identical  to 
itself  if  rotated  through  90°  about  its  axis.  Thus  the 
phase  space  emittance  of  the  particle  beam  in  the  focusing 
plane  is  congruent  with  that  in  the  defocusing  plane. 

From  this  it  follows  that  only  that  part  of  the  emittance 
common  to  both  the  ellipses  of  Fig.  3-1  is  accepted  and 
therefore  the  transverse  acceptance  at  the  plane  of  injec¬ 
tion  is  —  times  the  shaded  area. 

n 
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The  transverse  acceptance  at  any  other  plane  can 
also  be  found  as  follows: 

The  equations  in  matrix  form  between  the  plane  of  injec¬ 
tion  and  any  other  transverse  plane,  denoted  by  suffix  1, 
can  be  written  as 


for  vertical  motion 
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for  horizontal  motion 
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of  the  transfer  matrices  for  vertical  and  horizontal  planes 
respectively . 

Using  the  relations  a  vd  ^-b  -^c  ^  =  1  and  a.2<^2~^2c2  =  d 
the  foregoing  equations  become 
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(3-64) 


Using  the  boundary  values  of  (x0,x0)  and(yQ,y0)  given 
by  equations  3-59  and  3-60  the  above  equations  become 
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.  2  .  p 

( d px — b x -]  )  (a-.x-1-c-iX-i)^ 

p1  1  +  ■  1  1  l  1  •  =  1  - (3-65) 
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- -  +  - p -  =  i - (3-66) 

(-)  l—Y 
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By  rotation  of  axes  the  foregoing  equations  can  be 
written  as 

2  2 

— —  +  — —  =  1  - - - - - (3-67) 

A2b2  B2b2 

y2  +  y2  2  =  1  - (3-68) 

2  (b  \  2  2  fb\ 

c  (?)  D  hi 

where 

A=l//(d^  +  Ci32)COS24)x+(b^+ai3^)SIN^qx-2(d1b1  +  c1a132)SIN  4>XC0S  <J>X 
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and  where  the  axes  have  been  rotated  by  angles  cf>  and  <J> 

x  y 

given  by 


4> 

cl) 


X 


y 


“tan  1/2(b1d1+c1a1B2)/( (b2-d2 )+B2 (a2-c2) ) 

itan"1|2(b2d2+hi^_)/((b2-d|)+T-(a2-c|))'’ 


Equations  3-65  and  3-66  are  shown  in  Pig.  3-2  where 
the  acceptance  at  the  transverse  plane  denoted  by  sub¬ 
script^-  is  the  shaded  area  common  to  both  ellipses. 


Fig.  3-2  Transverse  Acceptance  at  any  Plane  Along 


the  Accelerator 
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CHAPTER  4 

NUMERICAL  STUDIES  OF  THE  MICROPARTICLE 

DYNAMICS  IN  LINEAR  PARTICLE  ACCELERATOR 

4.1  Introduction 

In  Chapter  2  of  this  thesis  the  basic  formulas  of 
the  particle  motion  in  the  accelerator  sections  and  in 
quadrupole  lenses  were  derived.  The  results  of  this 
study  were  applied  in  Chapter  3  to  investigate  the  pro¬ 
perties  of  particle  motion  in  the  transverse  planes  of  an 
idealized  accelerator  structure  in  which  all  sections  are 
identical.  The  study  of  Chapter  3  has  neglected  any 
velocity  gain  of  the  particles,  and  for  this  reason  the 
approach  that  has  been  used,  while  providing  useful 
physical  insight  cannot  be  employed  to  make  exact 
numerical  calculations . 

The  purpose  of  the  present  chapter  is  to  make  an 
exact  computer  study  of  the  dynamics  of  particles  in  an 
accelerator  structure.  The  formulas  used  here  are  those 
that  have  been  derived  in  Chapter  2. 


4.2  The  Computer  Program 


' 
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4.2.1  General  Discussion 

A  program  has  been  written  for  the  IBM  System/360 
computer.  The  purpose  of  the  program  is  to  design  the 
accelerator  drift  tube  structure  for  a  given  set  of 
nominal  particle  input  parameters „  These  nominal  input 
parameters  are  the  particle  charge  to  mass  ratio,  its 
injection  potential,  its  synchronous  phase,  the  operat¬ 
ing  frequency,  the  peak  voltage  in  the  gap,  the  radius 
of  the  drift  tubes,  and  the  ratio  of  the  gap  length  to 
the  drift  tube  length.  Once  the  structure  has  been 
designed,  its  properties  are  then  studies  while  the 
input  conditions  are  varied  about  their  nominal  values. 
This  study  for  varying  input  conditions  is  conducted 
first  without,  and  then  with  quadrupole  focusing  elements 
within  the  above  designed  structure.  The  results  are 
obtained  partially  in  tabular  form  and  partially  in 
graphical  form.  They  are  presented  throughout  the 
remainder  of  this  chapter, 

4.2.2  Function  of  the  Program 

The  dimensions  of  an  accelerator  structure  without 
quadrupoles  as  shown  in  fig.  4.1  are  computed  in  the  sub¬ 
routine  STRUCT  for  a  synchronous  particle  running  along 
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the  axis  of  the  accelerator.  This  subroutine  can 
also  be  used  for  designing  the  structure  for  an  off 
axis  synchronous  particle  although  this  has  not  been 
done  in  this  work.  Subroutine  STRUCT  also  computes 
the  synchronous  velocity  and  the  synchronous  energy 
at  the  various  drift  tube  centres. 


-1 


LSI ( 1)  -+ 
_  G(0) 

(w-  2 


* - LS  1(2) - » 


K-1— N 

i  I  i 


G  ( 1 ) 


k--^T 

•  I  i 


G  ( 2 ) 


i - 


'  “I  ■  r 


i  i 


—  L  D  ( 1 ) — —  LD(2)— *j 


LD  ( 3 ) 


LS  (1) 
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Injection  Plane 


Note:  The  same  symbols  are  used  m  the  program 


Fig.  4.1  The  First  Two  Sections  of  the  Accelerator 

The  main  program  is  divided  into  two  parts.  The 
first  part  is  for  the  structure  without  quadrupoles 
and  the  second  part  with  quadrupoles.  However,  the 
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second  part  can  also  be  used  for  the  structure  without 
quadrupoles  if  one  takes  the  quadrupole  voltage  to  be 
equal  to  zero. 

The  subroutine  RESULT  uses  the  dimensions 
found  by  STRUCT  to  compute  the  velocity  gain,  energy 
gain,  change  in  particle  phase  for  a  given  section 
for  nonsynchronous  particles .  Further  the  subroutine 
RESULT  computes  the  change  in  slope  of  each  particle 
trajectory  at  the  centre  of  the  gap  in  that  section. 
Using  the  findings  of  the  subroutine  RESULT,  the  main 
program  now  computes  the  actual  particle  trajectory. 
This  process  is  repeated  section  by  section  for  the 
complete  length  of  the  accelerator. 

The  subroutine  PLOTNG  produces  different  graphs 
from  the  digital  output  of  the  main  program.  A  list 
of  identification  of  the  variables  used  in  the  program 
is  given  below. 


Program  Variables 

Identification  of  the  Variables 

A 

Inner  radius  of  the  drift  tubes  in 

meters 

ALPHA (N) 

Same  as  k^  in  equation  2-4  of  page  8 

ANGL 

Angle  in  radians 

BI1 

Io(nr/Lsn> 

BI2 

I0(na/LSn) 
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Program  Variables 

Identification  of  the  Variables 

BI3 

I1(nr/LSn) 

BI4 

I2  (IIr/LSn) 

B15 

I1(na/LSn) 

BESI 

Library  subroutine  for  the  Bessel 

functions 

D 

Distance  of  the  quadrupole  tips 

from  the  axis  of  the  accelerator  in 

meters 

DEG 

The  projections  of  the  angle  of  the 

trajectory  at  injection  in  degrees 

on  both  horizontal  and  vertical  planes. 

DELTA 

Same  as  A^  in  equation  2-13  of  page  12 

for  synchronous  particle 

DELTAl 

Same  as  in  equation  2-13  of  page  12 

for  nonsynchronous  particle 

DELV 

Velocity  gain  in  Meters/Sec.  for  a 

nonsynchronous  particle  in  one  section 

DEL  VS 

Velocity  gain  in  Meters/Sec.  for  a 

synchronous  particle  in  one  section 

DENERG(N) 

Difference  between  the  synchronous 

and  the  nonsynchronous  energy  in  volts 

at  the  middle  of  the  nth  drift  tube 

DENG 

Energy  gain  in  volts  for  a  nonsynchronous 

particle  in  the  nth  section 

. 


- 
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Program  Variables 

Identification  of  the  Variables 

DEN  GY 

Energy  gain  in  volts  for  a  synchronous 

particle  in  the  nth  section 

DPHI (N+l) 

Difference  between  the  nonsynchronous 

and  the  synchronous  phase  in  degrees 

at  the  nth  gap  centre. 

DX,  DXPR , 

DXl 

Change  in  transverse  velocity  in 

meters/sec.  in  the  vertical  plane  and 

at  the  gap  centres 

DY ,  DYPR, 

DY1 

Change  in  transverse  velocity  in 

meters/sec.  in  the  horizontal  plane 

and  at  the  gap  centres 

EM 

Charge  to  mass  ratio  in  coulombs/Kg 

ENERGS (N) 

Energy  of  the  synchronous  particle 

in  volts  at  the  middle  of  the  nth 

drift  tube 

ENERGY (N) 

Energy  of  the  nonsynchronous  particle 

in  volts  at  the  middle  of  the  nth  drift 

tube 

F 

Operating  frequency  in  cycles/sec 

FA 

Ratio  of  the  axial  length  of  the 

quadrupoles  to  the  drift  tube  length 

FB 

Ratio  of  the  gap  length  to  the  drift 

tube  length 

G  (N) 

Length  of  the  nth  gap  in  meters 

. 


* 


. 
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Program  Variables 

Identification  of  the  Variables 

IK 

Dummy  variable 

IERI,  IER2,  IER3 , 
IER4,  IER5 

Error  Codes 

J2  (N) 

An  array  of  data  points  to  be  plotted 

for  graphs 

L 

Number  of  drift  tubes 

LSI ,  SS2 

Distance  in  meters  from  the  (n-l)th 

gap  centre  to  the  nth  gap  centre 

LD  (N) 

Length  of  the  nth  drift  tube  in  meters 

LS  (N) 

Length  of  the  nth  section  in  meters 

M 

A  variable  which  determines  the 

symbols  to  be  used  in  the  graphs 

NA 

Total  number  of  points  to  be  plotted 

in  the  graph 

PH 

Nonsynchronous  phase  in  degrees  at 

the  middle  of  the  first  drift  tube 

PHI  (1)  ,  PHIS  (1) 

Nonsynchronous  and  synchronous  phase 

in  radians  at  the  middle  of  the  1st 

drift  tube 

PHIS (N+l) 

Synchronous  phase  in  radians  at  the 

nth  gap  centre 

PHI (N+l) 

Nonsynchronous  phase  in  radians  at 

the  nth  gap  centre 

PHS ,  PHISl(N) 

Synchronous  phase  in  degrees  at  the 

nth  gap  centre 

» 
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Program  Variables 

Identification  of  the  Variables 

PI 

Same  as  k  defined  in  equation  2-25 

of  page  21 

Ql,  Qll,  Q12 

The  arguments  of  COS,  COSH,  SIN, 

SINH  of  equations  2-34  and  2-35  of 

page  23  and  24  respectively 

R(N) 

Radial  excursion  in  meters  at  the 

middle  of  the  nth  drift  tube 

Rl 

Radial  excursion  in  meters  at  the  gap 

centre 

RPRIME (N) 

Radial  velocity  of  the  particle  in 

meters/sec.  at  the  middle  of  the  nth 

drift  tube 

RPRIM6 

Average  radial  velocity  in  meters/sec. 

at  the  gap 

S,  SA,  S3,  S4,  S5 , 
S6  ,  TS,  TS1,  TS2 

The  times  taken  by  the  particle  in 

transversing  different  portions  of 

an  accelerator  section 

TOTLEN 

Total  length  of  the  structure  in  meters 

TT  (N) 

Transit  time  factor  for  the  particle 

at  synchronous  velocity 

VS  (N) 

Synchronous  velocity  in  meters/sec. 

at  the  middle  of  the  nth  drift  tube 

V  (N) 

Nonsynchronous  velocity  in  meters/sec. 

at  the  middle  of  the  nth  drift  tube 
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Program  Variables 

Identification  of  the  Variables 

VM 

Peak  voltage  in  the  gaps  in  volts 

VMl 

Voltage  of  the  particle  at  injection 

in  volts 

VM2 

Voltage  between  the  adjacent 

quadrupoles  in  volts 

V2  (N) 

Error  in  the  computation  of  the  length 

LSI (N) 

W 

Same  as  Q  in  equation  2-29  of  page  21 

X  (N) 

Transverse  excursion  in  meters  in 

the  vertical  plane  and  at  the  middle 

of  the  nth  drift  tube 

XA,  XI,  X2,  X3 ,  X4, 
X5  ,  X6 

Transverse  excursions  in  meters  in 

the  vertical  plane  and  at  different 

cross-sections  within  an  accelerator 

section 

XPRIME (N) 

Transverse  velocity  in  meters/sec. 

in  the  vertical  plane  and  at  the 

middle  of  the  nth  drift  tube 

XPRIM1,  XPRIM2 , 
XPRIM3,  XPRIM4, 
XPRIM5 ,  XPRIM6 

Transverse  velocity  in  meters/sec 

in  the  vertical  plane  and  at 

different  cross-sections  within  an 

accelerator  section 
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Program  Variables 

Identification  of  the  Variables 

Y(N) 

YA,  Y 1 ,  Y2 ,  Y3  f 

Y4 ,  Y5 ,  Y6 ,  YPRIMl , 
YPRIM2  ,  YPRIM3 
YPRIM4 ,  YPRIM5 , 
YPRIM6 

Same  as  X(N),  XA,  Xl ,  etc.  described 

above  bur  in  the  horizontal  plane 
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4.2.3  Flow  Charts 


MAIN  PROGRAM 


62 
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SUBROUTINE  STRUCT 


64 


SUBROUTINE  PLOTNG 
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4.2.4  Listing  of  the  Program 


COMMON  EM,VM1,F,R1,A, VM,L,FB,PH,DELTA1,LS(40) , 

1 G ( 40 ) , ENERGS ( 40 ) *  ENERGY ( 40 ) , DVEL <  4G ) , DPH I ( 40 ) , 
1VS ( 40 )  ,LD( 40 ) , PHI S (40 ) , PH  I (40 ) , ALPHA (40 ) , V( 40 ) , 
1 DENERG ( 40  ) 

DIMENSION  X ( 40) , Y(40) , XPRIME (40 ) , YPRIME( 40) , 
1XPRTM<40) , YPR I M ( 40 ) , R( 40 ) ,  J2  ( 40 ) »  DATA ( 1024) 

PFAL  LD»LS»J2 

CALL  PLOTS  ( DAT  A ( 1 )  ,4096) 

CALL  PLOT  (0,0,3.0,23) 

NOMINAL  INPUT  SPECIFICATIONS 


REAO( 5,40 )A, VM, VM1 , F , E M , L , D , F A , FB 
40  FORMAT (F5, 2  »3F9,1  ,F6,I  ,  I2,3F5.3) 


COMPUTATION  OP  STRUCTURE  DIMENSIONS 
CALL  STRUCT 

INJECTION  PARAMETERS  OF  NONS Y NC HR ONOU S  PARTICLE 
AT  THE  MIDDLE  OF  FIRST  DRIFT  TUBE 

R  E  AD ( 5 , 5  )  X (  l  ) ,  Y ( 1  ),DEG,PH 
5  FORMAT ( 4F10. 5 ) 

M=  I  30 

V(  I  )  =S QRT ( 2.*EM*VM1  ) 
XPRIME(1)=(TAN(0EG/57.29582))*V<1) 

YPR I  ME (1 )  =  ( TAN( DEG/57.295  82)  )*V( 1 ) 

FNERGY (1  )  =  .5*(  V( 1 ) 442 • ) /EM 

PHI  ( 1) =(-90.+PH  )/57.?9582 

DPHM  1  )  =  (  PHI  ( 1  )  -PHIS  ( 1  )  )  457.  2  958  2 

STRUCTURE  WITHOUT  QUADRUPOLES 

DO  80  N  =  1  ,  L 

COMPUTATION  OF  NONSY NCHRONOUS  PARAMETERS  AT 
THE  MIDDLE  QE  (N+1)TH.  DRIFT  TURE  FOR  STRUCTURE 
WITHOUT  QUADRUPOLES 

J2( N  )  =FL  OAT ( N) 

N  A  =  N 


FOP  THE  COMPUTATION  OF  PARTICLE  TRAJECTORY, 

EACH  SECTION  HAS  BEEN  DIVIDED  INTO  SEVFRAL 
SUB-SECTIONS.  THESE  SUB-SECTIONS  ARE  DEFINED 
IN  THE  FOLLOWING  COMMENTS  WHERE  THE  DESCRIPTION 
OF  S4,S5,TS,S,S3,TSi ,SA,S31,TS2  IS  ALSO  GIVEN 

S  4  IS  THE  TIMF  TAKEN  BY  THE  PARTICLE  TO 
TRAVERSE  THE  LENGTH  FROM  TH<=  CENTRE  OF  THE  NTH. 


, 
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DRIFT  TUBE  TO  THE  CENTRE  OF  THE  NTH.  GAP 

S4=.5*(1.+F3)*LD(N)/V(N) 

XA=X(N)+S4*XPRIME(N) 

Y A  =  Y ( N  )  +  S4AYPR I  ME ( N ) 

PI  =SQRT( ( XA**2  )+(YA**2) ) 

XPRIM1=XPRIME(N) 

YPR IM1=YPRI MF(N) 

CALL  RFSULTC  NtXPRIMl ,YPRIMlf XA, YAf DXPR»OYPR ) 
XPRIMF(N+1  )  =  XPR I  ME ( N ) +  DXPR 
YPR  I  me (N+l ) =YPR IME( N ) +OYPR 

S5  IS  THE  TIMF  TAKEN  BY  THE  PARTI CLF  TO  TRAVERSE 
THE  LENGTH  FROM  THE  CENTRE  OF  THE  NTH.  GAP  TO 
THE  CENTRE  OF  THE  IN+1ITH.  DRIFT  TUBE 

S5=( ( . 5*LD(N +1 )  )  +  < .5 *G(N) ) )/(  V(N+III 
X ( N+l )=XA+S5 *X PRIME (N  +  l ) 

Y(NU)=YA+S5*YPRIWE(N+1) 

XPRTM(N)=(ATAN(XPRIME(N)/V(N)))*57. 29582 
YPRIMIN )=(ATAN(YPPIME(N)/V(N)  )  ) *57.29582 
P(N+1)=SQRT((X(N+1)**2)+(Y(N+1)**2)) 

I F ( { R ( N  +  l  ) )  .GT.  A  )  GO  TO  7 
I F ( ABSIPHI (N+l ) )  .GT.  6.28318)  GO  TO  7 
8C  CONTINUE 

OUTPUT  IN  TABULAR  FORM 

7  WRITE  (6,25  )<(  K,  X(  K)  t  Y(  K  ).  DPHI  (  K  ),  DVEL  (K  >,  DENERGK  ))  ♦ 
1 K= 1 »  NA  ) 

25  FORMAT!  (5Xf  I?*5( 5X,E11.4)  )  ) 

CUTPUT  IN  GRAPHICAL  FORM 

CALL  P  LOT  NG  ( J 2 » N A  ,  X  ,  Y , L t M  ) 

STRUCTURE  WITH  AND  WITHOUT  QUADRUPLES 

COMPUTATION  OF  FIFLD  GRADIENT  FOR  QUAOPUPOLE 
VOLTAGE  VM2.  FOR  STRUCTURE  WITHOUT  QUADRUPOLES 
SET  V M2  =G • 

READ! 5,45 )VM2 
45  FORM AT (FIG. 2) 

P 1=  VM  ? / ( D**2.  ) 

W  =  SQR  T  (  P 1  <=  C  E^l  ) 

DO  8°  N  =  l »  L  ,  2 

COMPUTATION  OF  NON SYNCHRONOUS  PARAMETERS  AT  THE 
MIDDLC  OF  (N  +  l  ) T  H •  AND  (N+2)TH.  DRIFT  TUBES  FOR 
THE  STRUCTURE  WITH  AND  WITHOUT  QUADRUPOLES 

J  2 ( N) =FLOAT ( N) 

N  A  =  N 

TS  IS  THE  T I  ME  TAKEN  BY  THE  PARTICLE  TO  TRAVERSE 
THF  SECOND  HALF  LENGTH  OF  THE  QUADRUPOLES  IN 
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THE  NTH.  DRIFT  TUBE 

TS=(FA*LD(N) )/( V( N)*2. ) 

01=  W*T S 
02  =  COSH (01  I 
Q3=SINH(Q1 ) 

04=  COS (01  I 
05=S IN( 01 ) 

IR ( W  . EQ.  0. )  GO  TO  26 
Q6=Q5/W 
07=03/ W 
GO  TO  27 

THE  FOLLOWING  THREE  STATEMENTS  ARE  USED  FOR 
STRUCTURE  WITHOUT  QUADRUPLES 

?6  Q6=PA*LD(N) /(2.*V(N) > 

04  =  07 
D  =  A 

27  X  1=  (  X(N)*Q4)  •*•(  XPRIME  (N)*Q6) 

XPRIM1  =  -(X(N) *Q5*W)  4-  ( XPRI ME (N ) *Q4) 

Y1  =  {  Y  {  N  )  *02  )  4-  (  YPRI ME ( N) *07) 

YPRIM1  =  (  Y(N)*Q3*W)4-(  Y  PRIME  <N)  *02) 

S  IS  THE  TIME  TAKEN  BY  THE  PARTICLE  TO  TRAVERSE 
THE  LENGTH  FROM  THE  END  OF  THE  QUADRUPQLE S  IN 
THE  NTH.  DRIFT  TUBE  TO  THF  CENTRE  OF  THE  NTH.  GAP 

S=(G(N)+(( l.-FA >*LD(N)  )  )/(2.*V(N)  ) 

X2  =  X  1  -*-  (  S*  X  DR I  Ml  ) 

Y 2  =  Y 1  +  ( S^YDR  I  Ml  ) 

R1=S0RT (  ( ABS ( X2 ) ) **2  4- ( AB S ( Y2  )  )  **2  ) 

CALL  RESULT ( N, XPRIM1 f  YPRIM1 ,X2* Y2,DX ♦ DY ) 

XPR I M2=XPR I Ml+n X 
VPR I  M2  =  Y P R  I  Ml +0Y 

S^>  IS  THE  T I  MF  TAKEN  BY  THE  PARTICLE  TO  TRAVERSE 
THE  LENGTH  FROM  THF  CENTRE  OF  THE  NTH.  GAP  TO 
THE  BEGINNING  OF  THE  QUADPUPOLES  IN  THE  (N-H)TH. 
DRIFT  TUBE 

S3= ( G( N) +( ( 1 .-F A)*LD( N+l ) ) ) / ( 2. *V ( N+l ) ) 

X 3  =  X  2  +  (  S3*  XPR  I  m?  ) 

Y  3  =  Y  24-  (  S3*YPR  IM2  ) 

TS1  IS  THE  TIME  TAKEN  BY  THF  PARTICLE  TO  TRAVERSE 
THP  FIRST  HALF  LENGTH  OF  THE  QUADRUPOLFS  IN  THF 
(Nfl)TH.  DRIFT  TUBE.  THIS  IS  ALSO  THE  TIME  TAK<=N 
RY  THE  PARTICLE  TO  TRAVERSE  THE  SFCOND  HALF 
LENGTH  OF  the  (  N4-1  )TH.  QUADRUPOLFS 

T  S  1  =  (  (  FA*LD(  N4- 1  )  )  /  (  V  (  M4- 1  )  *  2  .  )  ) 

011=W*TS1 
021 =COSH ( 01 1 ) 

0?1=S INH( Oil ) 

0 4 1  =  C 0 S ( Oil  ) 

Q51=SIN(Q11 ) 
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I F ( W  .  EQ.  0. )  GO  TO  28 
Q6 1=Q  51/W 
Q7 1 =Q  3 1 / W 
GO  TO  2Q 

THE  FOLLOWING  TWO  STATEMENTS  ARE  USED  FOR 
STRUCTURE  WITHOUT  QUADRUPOLES 

28  Q61=FA*L0(N+1)/(2.*V(N+1)  ) 

Q7 1 =Q  61 

2 Q  X ( N+l )  =  ( X3*Q21 ) +  (  XPPIM2*Q71) 

XPR I  ME ( N  +  l )  =  ( X3*Q31*W)+(XPRIM2*Q21 ) 

XPR IM (N+l )  =  (  AT  AN( XPP I  ME (N  +  l ) /V(N+1 )  )  ) *5  7.  29  58? 
Y(N+1)=( Y3*Q41 )+( YPPIM2*Q61  ) 

P (N+l ) =S  QRT ( ( X( N+l )  *  *2 )  +  (  Y  {  N  + 1 >**2 ) ) 

I  F ( ( R ( N  +  l  )  )  .  GT.  D  )  GO  TO  15 

I F ( ABS( PHI ( N+l ))  .GT.  6.28318)  GO  TO  15 
YPRIME (N+l ) =-( Y3*Q51*W)+(YPRIM2*Q41 ) 

YPR IM (N  +  l )  =  (  AT  AN (  YPR  I  ME ( N  +  l ) / V ( N+ 1 )  )  ) *5  7.  29  58  2 
X4=X(N+1)*Q21 +XPRIME(N+1)*Q71 
XPPIM4=(X{N+1)*Q31*W)+(XPRIME(N+1)*Q21) 

Y4= ( Y (N+l ) *Q41 ) +( YPRIME (N+l ) *061 ) 

YPR I M  4=— (Y(N+1)*Q51*W)+(YPRIME(N+1)*Q41) 

SA  IS  THE  TIME  TAKEN  BY  THE  PARTICLE  TO  TRAVERSE 
THE  DISTANCE  FROM  THE  END  OF  THE  QUADRUPOLES  IN 
THE  (N+l ) TH .  DRIFT  TUBE  TO  THE  CENTRE  OF  THE 
(N+l ) TH.  GAP 

S  A  =  (  G  (  N  + 1 )  +  (l.-FA)*LD(N+l)  )/(2.*V(N+l)  ) 

X  5  =  X4  +  (  SA*XPRIM4 ) 

Y  5  =  Y  4+ (  SA*YPP I M4 ) 

P  1  =  SQRT(  ( ARS ( X5  )  )**2  + ( AB S ( Y 5 )  ) ** 2  ) 

CALL  RESULT ( N  +  l  , X PR  I M4 , YpP I M4 , X 5 , Y 5 , DX 1 , DY1  ) 

XPR  IM5=XPRI M4+DX1 
YPR  I  m 5  =  Y P R I M4+D Y1 

S 3 1  IS  THE  TIME  TAKEN  BY  THE  PARTICLE  TO  TRAVERSE 
THE  LENGTH  FROM  THE  CENTRE  OF  THE  (N+1ITH.  GAP 
TO  THF  BEGINNING  OF  THE  QUADRUPOLES  IN  THE 
( N  +  2  )  T H .  DRIFT  TUBE 

S31=(G(N+1)  +  (1.-FA)*LD(N  +  ?)  )/(2.*V(N+2)  ) 

X6  =  X5  +  S31*XPRTMc; 

Y6=Y5+S31 *YPRI M5 

TS2  IS  THE  TIME  TAKEN  BY  THE  PARTICLE  TO  TRAVERSE 
THE  FIRST  HALF  IENGTH  OF  THE  QUADRUPOLES  IN 
THF  ( N+2 ) TH.  DRIFT  TUBE 

TS?  =  ( FA*LD<  N+? )  )/ ( V ( N  +  2 ) *2.  ) 

Q 1 2  =  W*  T  S? 

Q??=COSH( Q12 ) 

Q?2=SINH(Q1?) 

Q42-C0S ( Ql? ) 

Q52=SIN(Q1? ) 

I F ( W  .EQ.  0.)  GO  TO  50 


69 


Q62=Q52/W 
Q7  2  =  0  32/ W 
GG  TG  54 

THE  FOLLOWING  TWO  STATEMENTS  ARE  USED  FOR 
STRUCTURE  WITHOUT  QUADRUPOLES 

5G  Q62  =  FA*LD(  N  +  2) /(  2.*V<  N  +  2)  ) 

Q7  2  =  062 

54  X(N+2)=(X6*Q42)+(XPRIM5*Q62) 

XPRIMEIN  +  2  )=-(  X6*Q52*W )♦(  XPRIM5*Q42) 

XPRIM(N+2)=(  AT  AN! XPR IMF (N+21/VCN+2 ) )  )*57. 29582 
Y(N+2)=(Y 6*022 ) +< YPR I M5* 07 2 ) 

YPRTME(N+2)=(Y6*Q32*W  )  4-  (  YPR  I  M5*Q22  ) 

YPRIM<N+2) =( AT AN( YPR I ME (N+2 )/V(N+2 )) )*57. 29582 
N  A  =  N  +  i 

J  2(  N+l  )  =  F  L  0  A  T  (  N  + 1  ) 

P (N+2) =SQRT< ( X( N+2 ) ^  ^2 ) +( Y(N+2 ) **2 ) ) 

IF ( ( P (N+2 ) )  • GT •  0  )  GO  TO  15 

I F  M  ABS ( PH  I ( N+2 ) ) )  .GT.  6.28318)  GO  TO  15 
89  CONTINUE 

OUTPUT  IN  TABULAR  FORM 

15  WRITE! 6, 91 ) (  (K, X(K  ) , Y(K ) ,DPHI ( K) ,DVEL  <K ) ,DENERG(K  > )  , 
1 K= 1 ♦ NA ) 

Q 1  FORM  AT!  ( 5X , I  2, 5! 5X, El  1.4) ) ) 

OUTPUT  IN  GRAPHICAL  FORM 

CALL  PLOTNG  ( J 2 » N A , X , Y , L , M  ) 

CALL  PLOT  (0.0,0.0*099) 

STQR 

END 
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3HNI1N03  29 

((9d  +  *Z)/HN)01j}t9d))-((8=)^9,+  *I)/2SS)  =  (T  +  N)Q_l 

2  S  S  =  (  I -»-N  )  I  S 1 
1 9  Gi  03  (2.0-3  *  T  *1V  TA  )  31 

(  I +  N  )  I  S”l-2 S  S  =T  A 
(  3*  *  2 )  / (  T+N ) SA=2SS 
SA"13G+(N  )SA=(  I+N  )SA 
( N  )  11* (  (N)SA/(  (T+N)SIHd)S03*WA*W3  )=SA“I3G 

2I8/I ia*^a*(N)ii 
( 2831 *  2 1  8  *  0  *  2C ) I S38  T1V3 
(  1831  4 118*0* 10  )  IS38  IIVD 

ea/ua  >ni  s=*?g 

(  IN) SI* *2)/ (  ( N  )0*b9  I*?  I *t  )  =  £0 
( N  >  $  1/  V*o&IVr£  =  2Q 
(N)S1/T8*69lVT,£=ia 
(N)3+(  (  (I+N)Cn+(N)a“l>*S#)*(N)Sl 

02*x=r  29  ua 

S3  OV Id 

1VWI330  1  Gidn  1038803  (I+N)IS1  GNV  (I+N)(]l  30 
S  3  H  *1 V  A  3H1  S3indW03  NO  I IV  831 1  3NIMGT1G3  3H1 

(N)  IS~I=(I  +  N)TS1 
( N  )CJO=  ( I+N  )G1 
(N  )  3W  I  dd 8 v9S  +  (  N)  8  =  18 
(  N  )  S A /  C  (83+M  )*(N)ai*9#  )  =  9S 
(2  )SJHd=(l  +  N) SIHd 
(N  )0~U83=  (  N  )  0 

"l  *  I  =N  oi  oa 

28962 r19/SHd= ( 2)SlHd 
28962*19/( SHd+#06-)=( T) SIHd 
(  (  9  3+  *  I  )/(  I  )  ISD  =  (T  )  Q  1 
(3**2)/( (  I)SA)  =  (I)IS1 
*0=N31101 

W3/ (  #2**( I  )SA)*9 *  =  ( T  )  S383N3 
( IwA*W3**2  ) 180S= ( I  )  S  A 
•0=(  I  )  3 W i 8d  8 

•o=n )  8 

(  1*93  11VW8  03  2  *7 
SHd(2*7*9)GV38 

svavisi  "IV 3 8 

( 0 V ) 3w I 8d  81 

*  (0*7  )  8*  (0*7  )li*  (0*7  )2A  *  (  0*7)  IS  IHd*  (0*?)  IS!  NO  I SN  3W  IG 

(  0*7)083N3GI 

*(0*7)A*(0*7)VHdlV*(0*7)IHd*(0*7)SIHd*(0*7)Ql*(0*7)SAI 

*  (0*7)  IHda*  (0*7)  13Aa*  (0*7)A083N3*  (  0*7)S083N3*  (  0*7  )0I 

*  (0*7)  S"l*TV1130*Hd  *83  *1*HA*V‘T8*3*  I  to  A  *  w  3  NO  In  WO  3 

1 3 0 8 1 S  3Nim088nS 

SNQISN3W10  3  8ni3n8lS  3H1  S3iridW03  SNIinOdSOS  SIH1 
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53  V2(N)=V1 

LS1!N+1 )=SS2 

LD(N  +  l»=(SS2/(l.+.5*FBn-(  !FB*LD(N)  )/!2.+FB>  ) 

IS  (N)  =  !.5*!LD(N)-*-LD(N+l)  )  )+G(N) 
TOTLEN=TOTLEN+LS(N) 

I  F  <  R  1  .EG.  0.  )  GO  TO  92 
CALL  BESI (D1 ,1 ,BI3,  IER3 ) 

C 2=D4*B  I3/B 12 

DELTA=(  E”*VM*S  IN(PHIS(N+1M/VS(N)  )*<C2/R1) 

GO  TO  93 
92  DELTA=0. 

9  3  S4  =  .  5*<  LD(  N  +  l  )+G(NI»/IVS{N+in 

R(N+1>=PH-(S4*!-DELTA*R1  +  RPRIME(N)  )  ) 
RPRIME(N+1)=(-DELTA*R1 )+«PRIME(N) 

OENGY=VM*(  (CGS<PHIS(N+1)  )*TT<NM-(RPRIME(NI* 
lSIN(PHIS<N-H))*C2/VS(Nm 
ENERGS ( N+l I=ENERGS(N )+DENGY 
PHTS1!N1=PHIS(N)*57.29582 
70  CONTINUE 

WR I TE ( 6t 30  )  (  ( K,LD( K ) ,G(K ) ,TT( K)  , VS ( K)  , ENERGS ( K ) ) , 
1 K=1 tL ) 

30  FORMAT! <5X, 12,5! 5X, Ell. 4)  )) 

WR ITE( 6,78 )TOTLEN 
78  FORMAT! /5X , Ell. 4) 

RETURN 

END 
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THIS  SUBROUTINE  COMPUTES  NONSYNCHRONOUS  VELOCITY, 
PHASE, ENERGY  AND  THE  CHANGE  IN  SLOPE  OF  THE 
TRAJECTORY  IN  THE  GAP 

SUBROUTINE  R ESULT ( N, X  PR  I  Ml , Y PR  I  Ml , X A , YA , DXPR , DYPR ) 
COMMON  EM,VM1,F,R1,A,VM,L,FB,PH,DELTA1,LS(4C)  , 

1 G (  40 ) , ENERGS (40 ) »  FNFRGY( 40 ) , DVEL (40) ,DPHI (40) , 
1VS(40 ) , LD( 40 ) , PHIS (40 ) , PHI (40) , AL PH A ( 40 ) ,V( 40) , 
1DENERG( 40 ) 

DIMENSION  TT ( 40 ) 

PEAL  L S  ,  L D 

ALPHA(N) =( VS ( N ) /V ( N ) )-l. 

IF  ( N  .  GT.  1 )  GO  TO  81 

PHI  (2)=( (PH) / 57. 29582  )  +  (  (  90. /5  7.  2958  2  )  *  A  L  PH  A  (  1  )  ) 

GO  TO  82 

81  PHI (N+l )=PHI (N) +3 .14159* ALPHA* N) 

82  DPH  I  (  N 1 )  =  (  PH  I  (  N+l ) -PH  I S ( N+l ))  *5  7.  2  9582 
Dl  =  3. 14159*R1/LS(N) 

D2=3. 141 59*A  /LS(N) 

D3= ( 3. 14159*G( N)  )/<2.*LS(N)) 

D4=S I N( D3)/D3 

CALL  BESI (D1 ,0,BI1, IERI ) 

CALL  8ESI(D2»0,BI2» IFR2) 

CALL  BESI ( D1,1,BI3* IER3 ) 

CALL  BESI  (01 ,2 ,BI4,  IER4 ) 

CALL  BESI ( D2 , 1 , B I  5 ,  IER5 ) 

D5=( ((( .7853975 /LS (N ) )*(G(N) **2. )*C0S(D3) ) 
1-(.5*G(N)*SIN(D3) ) )  /<D3**2. ) ) 

TT( N )=D4*BI 1/BI2 
C  2  =  D4*  B  1 3 / B  I  2 

C3=-(3.14159/(LS(N)**2. ))*(*( D4*(  R 1  *  B I 2*B I  3 
1-A*BI 1*BI 5  )  ) /( BI 2**2 • ))  +  ( ( BI 1/BI  2) *0  5 ) ) 

C4=- ( 3.14159/(LS(N)**2.))*((D4*((((BI2*LS(N) 

1/3. 14159)*( BI3+D1*B 14 ) )-A*BI3*BI5) /( B  I  2 **2 • 

1 ) ) )+( ( BI3/B I2)*05) ) 

A3=TT  (  N)  4-  AL  PHA  (  N)  *LS(N)*C3 
A5=C2-ALPHA(N)*LS(N) *C4 
DELV=(EM*VM*COS(PHI (N+l ) )/V(N) )*A3 
V( N+l ) =V( N ) +DEL V 
P VFL (N)=V(N)-VS(N) 

I F ( R 1  .LT.  l.E-07)  GO  TO  14 

DELTA1=( FM*VM*SIN(PHI ( N+l ) ) /V(N) ) * ( A5/R1 ) 

COMPUTATION  OF  THE  CHANGE  IN  TRANSVERSE  VELOCITIES 
IN  THE  VERTICAL  AND  THE  HORIZONTAL  PLANES 
DENOTED  BY  DXPR  AND  DYPR  RESPECTIVELY 

91  I F ( (  ( ABS ( XA  )  )  .GT.  l.E-07)  . AND. (  ( ABS ( Y A )  ) 

l.GT.  l.E-07) )  GO  TO  11 
I F (  (  ( AB  S ( XA  )  )  .GT.  l.E-07)  . AND.  (  ( ABS ( Y A )  ) 
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l.LT,  1.E-Q7) )  GO  TO  12 

GO  TO  13 

11  ANGL=ATAN(  ( ARS ( YA)  )/  (ABS(XA))) 

GO  TO  15 

12  DXPR=-(DELTA1*R1 ) 

D  Y  P  R  =  0  •  0 

GO  TO  7 

13  DXPR=0* 

DYPR--(DELTA1*R1  ) 

GO  TO  7 

15  DXPR--(DELTA1*R1*C0S{ANGL)  ) 
DYPR=-(DELTA1*R1*SIN(ANGL) ) 

7  IF(XA  .LT.  O.O)  GO  TO  16 
GO  TO  17 

16  DXPR=-DXPR 

17  I F ( YA  .LT.  0.0)  GO  TO  18 
GO  TO  IP 

18  DYPP=-DYPR 
GO  TO  I® 

14  DEL  T  A1  =  0. 

0XPR=0. 

DYPP  =0. 

RPRIM6=SQRT(  (XPRIM1**2  )  +  <YPRIMl**2)  ) 

GO  TO  20 

19  XPRI M6-XPRI Ml+DXPR/2. 

YPRIM6=YPRIMl+DYPP/2. 

RPR  I  M 6>  —  (  (  XA*X  PR  IM6  )  /  R1  )  +  (  (  Y  A*Y  PR  I  M6  )  /  R  1  ) 

20  DENG=VM*( <COS( PHI (N+l )  )*A3)-((RPPTM6/V(N)  >* 
IS  IN ( PHI ( N  +  l )  )*A5)  ) 

ENERGY ( N  +  l  )  =  ENF RGY ( N  )  +  DENG 
OENERGt  N ) = ENERGY! N)-ENERGS(N) 

RETURN 

END 
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THIS  SUBROUTINE  PLOTS  GRAPHS  FOR  TRANSVERSE 
EXCURSIONS  IN  THE  VERTICAL  AND  THE  HORIZONTAL 
PLANES. 


SUBROUTINE  PLOT NG ( J2 , N A , X  , Y , L  ,  M ) 

DIMENSION  J2  (40)  , X(40 ) , Y(40 ) 

REAL  J2 
DC  =5. 

CD  =  4. 

CPC  =CD  +  1 . 

J2  IS  THE  ABSCISSA  ARRAY  NAME  DENOTING  THE 
MIDDLE  OF  EACH  DRIFT  TUBE  AND  X,Y  ARE  THE 
ORDINATE  ARRAY  NAMES  DENOTING  EXCURSIONS  IN 
THE  VERTICAL  AND  THE  HORIZONTAL  PLANES 
RESPECTIVELY 

SET  UP  SCALE  FOR  X  ,Y  AND  J2 
J2 ( NA  +  1 ) =1  . 

J2(NA*2)=(FLQAT(L ) ) /DC 

SCI =-.03 

SC  2  =  • 06 

X  (  N  A  + 1  )  =  S  C 1 

X  (  N A  +  2 )  =  ( S C2 /CD  ) 

Y (NA+1 ) =SC1 
Y ( N A+2 ) = (  SC  2/C.D  ) 

DRAW  THE  AXES  FOR  THE  VFPTICAL  EXCURSIONS 

CALL  AXIS(0. 0,0.0,'  ' , -1 , DC , 0. 0 , J 2 ( NA+1 ) , 

1 J2(NA+2) ,20.0) 

CALL  AX  IS (0.0, 0.0,  '  '  , 1 , C D , 90 . 0 , X ( N A +  1 )  , 

1 X ( NA+  2 )  ,20,0) 

DRAW  THE  TRAJECTORY  FOR  VERTICAL  EXCURSIONS 
CALL  L I NF ( J  2 , X , NA , 1 ,1 ,M) 

THIS  SETS  THE  ORIGIN  FOR  THE  NEXT  GRAPH 
CALI  PLOT (0. 0,CDC ,-3 ) 

DRAW  THE  AXES  FOR  HORIZONTAL  EXCURSIONS 

CALL  AXIS(0. 0,0.0,'  •  ,-l,DC,0.0, J2(NA  +  1  )  , 

1 J  2 ( N A+2), 20.0) 

CALL  AXIS (0.0, 0.0,  '  • , 1 , CD , 9C . 0 , Y ( NA ♦  1  ) , 

I  Y ( NA  +  2 ) , 20 . 0 ) 


DRAW  THF  TRAJECTORY  FOP  HORIZONTAL  EXCURSIONS 

CALI  L INE ( J2 ,Y, NA, 1 , I , M) 

CALL  PLOT (0. 0,-CDC,-3 ) 

RETURN 

END 
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4.3  Computer  Output  and  Discussion  of  Results 

4.3.1  Structure  Dimensions  and  the  Synchronous  Parameters 

The  accelerator  structure  dimensions,  transit 
time  factor,  velocity  and  the  energy  of  the  synchronous 
particle  for  each  section  have  been  computed  for  the 
following  nominal  input  specifications.  The  result  is 
listed  in  table  1. 

Input  specifications: 

Inner  radius  of  the  drift  tube  =  4  cm 
Peak  operating  voltage  in  the  gap  =  150  Kv 
Voltage  of  the  particle  at  injection  =  300  Kv 
Operating  frequency  of  the  structure  =  30  KHz 
Charge  to  mass  ratio  of  the  particle  =  30  coulombs/Kg 
Synchronous  phase  at  the  gap  centre  =  -30° 

Ratio  of  the  gap  length  to  the  drift  tube  length  =  .25 
The  synchronous  particle  travels  along  the  axis  of  the 
accelerator . 

4.3.2  Particle  Trajectories 


Fig.  4.2  shows  the  trajectories  of  two  nonsynchronous 
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Structure  Dimensions  and  Synchronous  Parameters 


N  Length  of  Nth , 

Length  of 

Transit  time 

Synchronous 

Synchronous 

drift  tube 

Nth.  gap 

factor 

vel oci ty 

energy 

(cm) 

(cm) 

(Km/Sec) 

(Mv) 

1 

5.657 

1.414 

.5317 

4.243 

.3 

2 

6.380 

1.595 

.5929 

4.731 

.3691 

3 

7.024 

1.756 

.6407 

5.219 

.4461 

4 

7.661 

1.915 

.6806 

5.698 

.5293 

5 

8.280 

2.070 

.7136 

6.163 

.6177 

6 

8.879 

2.220 

.7411 

6.615 

.7104 

7 

9.460 

2.365 

.7642 

7.051 

.8067 

8 

10.020 

2.505 

.7837 

7.474 

.9060 

9 

10.560 

2.641 

.8004 

7.882 

1.0080 

10 

11.090 

2,772 

.8147 

8.278 

1.1120 

11 

11.600 

2.900 

.8271 

.  8.662 

1.2180 

12 

12.090 

3.024 

.8379 

9.034 

1.3250 

13 

12.570 

3.144 

.8474 

9,395 

1.4340 

14 

13.040 

3.261 

.8559 

9.747 

1.5440 

15 

13.500 

3.374 

.8634 

10.090 

1.6550 

16 

13.940 

3.485 

.8701 

10.420 

1.7670 

17 

14.370 

3.593 

.8762 

10.750 

1.8800 

18 

14.800 

3.699 

.8816 

11.070 

1.9940 

19 

15.210 

3.802 

.8866 

11.380 

2.1090 

20 

15.610 

3.903 

.8911 

11.680 

2.2240 

21 

16.010 

4.002 

.8953 

11.940 

2.3400 

22 

16.400 

4.099 

.8991 

12.270 

2.4560 

23 

16.780 

4,190 

.9025 

12.550 

2.5730 

24 

17.150 

4.287 

.9058 

12.830 

2.6900 

25 

17.520 

4.379 

.9088 

13.110 

2,8080 

26 

17.870 

4,469 

.9115 

13.380 

2.9260 

27 

18.230 

4.557 

.9141 

13.640 

3.0440 

28 

18.580 

4.644 

.9165 

13.910 

3.1630 

29 

18.920 

4,729 

.9188 

14.160 

3.2820 

30 

19.250 

4,814 

.9209 

14.420 

3.4010 

Total 

length  of  the 

structure  = 

:  5.05  Meters 

Note: 

Input  specification  for 

the  design  of 

the  structure 

is  given 

in  section  4 

.3.1. 

• 

‘ 
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particles,  with  same  charge  to  mass  ratio  and  injection 
voltage  as  the  synchronous  particle,  in  the  drift  tube 
structure  designed  above ,  The  particles  are  injected 
at  the  middle  of  the  first  drift  tube.  One  particle  is 
injected  on  the  axis  with  a  slope  such  that  the  projections 
of  the  initial  trajectory  on  the  horizontal  and  the 
vertical  planes  make  angles  of  0.2  degrees  with  the  axis 
of  the  accelerator.  The  other  particle  is  injected 
paraxially  with  an  initial  horizontal  and  vertical  dis¬ 
placement  of  .2  cm  each.  It  is  seen  that,  in  the  structure 
without  quadrupoles ,  the  particle  injected  on  the  axis 
collides  with  the  18th  drift  tube  wall  and  the  particle 
injected  paxaxially  collides  with  the  14th  drift  tube 
wall.  On  the  other  hand,  if  quadrupoles  are  used  both 
the  particles  traverse  the  complete  structure  of  30  drift 
tubes  without  being  intercepted. 

It  can  also  be  seen  from  fig.  4.2  that  the  tra¬ 
jectories  of  the  particle  for  a  structure  with  quadrupoles 
have  small  oscillations  with  a  period  of  two  section 
length.  The  amplitudes  of  these  oscillations  increase 
with  the  increase  in  transverse  excursion  of  the  particle. 
The  reason  for  this  behaviour  is  that  the  individual 
quadrupole  action  is  much  stronger  than  the  overall 
focusing  effect  and  that  as  the  excursions  increase  the 
focusing  and  the  defocusing  forces  become  more  intense 
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Note:  N  is  the  middle  of  the  Nth.  drift  tube. 

The  projections  of  the  slope  of  the  trajectory  A  at 
injection  in  both  horizontal  and  vertical  planes  are 
.2  degree  each.  The  particle  for  trajectory  x  is  in¬ 
jected  paraxially.  All  other  injection  conditions  are 
synchronous . 


FIG. 4. 2  Transverse  Excursions  With  and 

Without  Quadruples 
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(refer  to  equations  2-26  and  2-27) 

4.3.3  Phase  and  Energy  Oscillations 

In  fig.  4.3  the  phase  oscillations  and  the 
energy  oscillations  for  particles  whose  injection  phases 
differ  from  the  synchronous  phase  have  been  plotted. 

The  structure  used  here  is  the  one  designed  in  section 
4.3.1.  The  particle  is  assumed  to  move  along  the  axis 
of  the  accelerator. 

Fig.  4.3  shows  that  the  phase  slips  in  one  direc¬ 
tion  or  the  other  i.e.,  the  slope  of  the  phase  trajectory 
is  positive  or  negative  depending  on  whether  the  energy 
is  above  or  below  the  synchronous  energy.  It  is  noticed 
that  when  the  phase  is  synchronous  the  slope  of  the 
energy  (difference  between  nonsynchronous  and  the 
synchronous  energy)  trajectory  is  zero  implying  that  the 
gain  in  energy  of  the  nonsynchronous  and  the  synchronous 
particle  is  the  same. 

Fig.  4.4  shows  the  effect  of  different  injection 
energies  on  the  phase  oscillations .  It  is  seen  that  the 
damping  of  phase  oscillation  becomes  slower  as  the  injec¬ 
tion  energy  increases  above  the  synchronous  injection 
energy.  This  is  because  the  velocity  change  for  a  given 
change  in  energy  decreases  as  the  particle  energy 
increases.  Hence  A 4)  and  thus  the  rate  of  phase  damping 
decreases  as  the  energy  increases. 
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Note:  N  is  the  middle  of  the  Nth.  drift  tube.  The  ordinate 
above  or  below  the  horizontal  lines  through  the  zero 
phase  and  zero  energy  indicates  that  the  nonsynchronous 
phase  and  energy  are  greater  than  the  synchronous  phase 
and  energy  by  that  amount.  The  injection  parameters, 
except  for  the  injection  phase,  of  the  particles  are 
the  same  as  those  of  the  synchronous  particle.  Trajec¬ 
tories  with  the  same  symbol  in  phase  and  energy  plots 
correspond  to  the  same  set  of  injection  conditions. 


FIG.  4.3  Phase  and  Energy  Oscillations  for 
Different  Injection  Phases  and 
Same  Injection  Energy. 


81 


O 

O 


o 

o 

o 


Note:  The  injection  parameters,  except  for  the  injection 
energy,  are  the  same  as  those  of  the  synchronous 
parti  cl e . 


FIG.  4.4  Phase  and  Energy  Oscillations  for 
Different  Injection  Energies  and 
Same  Injection  Phase. 
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4.3.4  Range  of  Injection  Angle  and  Injection  Phase  for 
Accepted  Particles 

The  shaded  area  of  fig.  4.5  shows  the  region  of 
accepted  injection  phase  and  the  corresponding  accepted 
injection  angle  at  (0°-90°)  synchronous  phase.  Similarly 
the  boundaries  of  the  acceptance  parameters  for  the  syn¬ 
chronous  phases  (-15°-90°) ,  (-30°-90°)  and  (-45°-90°)  are 

also  shown  (note  that  the  factor  -90°  is  included  because 
the  particle  is  injected  at  the  middle  of  a  drift  tube) . 

It  is  clear  from  the  fig.  4.5  that  the  allowable  injec¬ 
tion  angles  have  been  greatly  improved  for  the  structure 
with  quadrupoles  over  the  structure  without  quadrupoles. 

It  is  also  seen  that  the  decrease  of  synchronous  phases 
from  (-0°-90°)  up  to  (-45°-90°)  is  accompanied  by  a 
decrease  in  acceptable  injection  angles  corresponding  to 
these  particular  injection  synchronous  phases.  This 
result  is  due  to  increased  radial  defocusing  forces  (refer 
to  equation  2-8)  in  the  gaps.  On  the  other  hand,  if  par¬ 
ticles  are  injected  with  an  angle  of  zero  degree,  then 
the  range  of  acceptable  injection  phases  increases  as  the 
synchronous  phase  is  decreased.  This  is  due  to  the 
increased  axial  focusing  forces  (given  by  equation  2-23) 
and  due  to  the  absence  of  a  radial  defocusing  force. 

It  may  also  be  noticed  that  the  acceptable  injec¬ 
tion  angle  gradually  increases  to  a  peak  value  as  the 
injection  phase  is  increased  from  the  injection  synchronous 
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Note:  The  particle  is  injected  on  the  axis  of  the  accelerator. 
$  is  the  synchronous  phase  at  the  injection  point  i.e. 
al  the  middle  of  the  first  drift  tube.  The  injection 
angle  is  the  angle  of  the  trajectory  of  the  particle  at 
the  point  of  injection  projected  in  the  horizontal  and 
the  vertical  planes. 

FIG.  4.5  TheRange  of  Injection  Phase  and 
Injection  Angle  at  the  Axis  for 
Accepted  Microparticles. 
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phase  and  then  very  rapidly  drops  down  to  zero.  This 
occurs  because  the  radial  defocusing  force  decreases 
as  the  injection  phase  is  increased  from  the  injection 
synchronous  phase.  When  the  injection  phase  crosses 
the  (-0°-90°)  injection  phase  the  net  radial  focusing 
force  in  the  gap  becomes  positive  and  hence  the  angle  of 
injection  increases.  However  at  the  same  time  the  axial 
focusing  effect  becomes  weaker,  and  eventually  as  the 
injection  phase  crosses  (-90°)  ,  the  axial  field  becomes 
defocusing.  Thus  if  the  injection  phase  is  further 
increased  the  phase  very  quickly  slips  by  greater  than 
360°  and  the  injection  angle  falls  to  zero.  It  may 
also  be  seen  from  fig,  4.5  that  as  one  decreases  the 
injection  phase  below  the  injection  synchronous  phase 
the  radial  defocusing  force  at  the  gap  increases  and 
thus  the  injection  angle  decreases.  As  the  injection 
phase  is  further  decreased,  the  axial  force  causes  the 
particle  to  execute  larger  and  larger  oscillations  about 
the  synchronous  phase.  Consequently  the  particle 
alternates  in  phase  such  that  the  radial  defocusing  force 
is  alternately  very  much  stronger  and  then  very  much 
weaker  than  the  radial  defocusing  force  at  the  syn¬ 
chronous  phase.  It  is  believed  that  the  net  effect  of 
this  large  amplitude  of  oscillations  leads  to  an 
increased  acceptance  of  injection  angles.  However  if  the 
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injection  phase  is  further  decreased  the  initial 
defocusing  forces  become  so  large  that  the  particle  will 
immediately  strike  the  structure  wall  and  thus  the 
injection  angle  falls  to  zero. 

4.3.5  Effect  of  Variable  Charge  to  Mass  Ratio  on  Particle 

Trajectories  for  the  Structure  Designed  in  Sec.  4.3.1 

It  is  seen  from  fig.  4.6  that  in  the  structure 
without  quadrupoles  the  particle  experiences  more  and 
more  defocusing  force  as  the  charge  to  mass  ratio  increases 
beyond  the  design  value  of  30  coulombs/Kg „  This  may  also 
be  seen  from  equation  2-8  which  shows  that  the  defocus¬ 
ing  force  in  the  gaps  increases  with  an  increase  in  charge 
to  mass  ratio. 

From  fig.  4.7  it  is  seen  that  in  the  structure  with 
quadrupoles  the  overall  focusing  force  in  the  first  few 
sections  diminishes  with  increase  in  charge  to  mass  ratio. 
But  after  these  first  few  sections  the  overall  focus¬ 
ing  force  increases  rapidly  with  increase  in  charge  to 
mass  ratio.  This  is  due  to  the  fact  that  increase  in 
charge  to  mass  ratio  increases  the  focusing  force  of  the 
quadrupoles  (refer  to  equation  2-28,  2-29).  In  the  first 
few  sections  the  increase  in  defocusing  force  in  the  gap 
is  larger  than  the  increase  in  focusing  force  of  the 
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Note:  The  particle  is  injected  on  the  axis  at  an  angle  of 
.1  degree  in  the  horizontal  and  the  vertical  planes. 
The  injection  phase  is  synchronous. 


FIG.  4.6  Transverse  Excursions  of  Particles 
With  Different  Charge  to  Mass 
Rations  for  Structure  Without 
Quadruples  Nominally  Designed 

for  JL  =  30  coulombs/Kg. 
m 


EXCURSION  IN  VERTICAL  PLANE  EXCURSION  IN  HORIZONTAL  PLANE 

(  In  Meter  x  100  )  (  |n  Meter  x  100  ) 

-0.30  -0.10  0.0  0.10  0.30  0.50  -0.30  -0.10  o.o  0.10  0.30  0.50 
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Note:  Same  note  as  in  Fig.  4.C 

FIG.  4.7  Transverse  Excursions  of  Particles 
With  Different  Charge  to  Mass 
Ratios  for  Structure  with  Quadrupoles 
Nominally  Designed  for 

A.  =  30  coulombs/Kg. 
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quadrupoles .  But  at  the  later  sections,  because  of  the 
increase  in  velocity,  the  increase  in  defocusing  force 
in  the  gaps  is  very,  very  small  compared  to  the  increase 
in  focusing  force  of  the  quadrupoles. 

It  has  been  seen  that  for  a  particle  moving  along 
the  axis  of  the  accelerator  the  range  of  acceptable 
charge  to  mass  ratio  is  26.2  to  83.5  coulombs/Kg.  If 
the  particle  is  injected  at  a  slope  of  .2  degree  each 
with  the  horizontal  and  the  vertical  planes  the  range  of 
acceptable  charge  to  mass  ratio  is  also  26.2  to  83.5 
coulombs/Kg  in  a  structure  with  quadrupole  voltage  of 
40  Kv.  But  for  the  structure  without  quadrupoles  the 
range  of  charge  to  mass  ratio  is  zero  for  a  particle 
injected  at  .2  degree. 

4.3.6  Effect  of  operating  Frequency  on  Structure  Design 

Additional  tables  of  accelerator  structure  dimens¬ 
ions  have  been  computed  for  an  operating  frequency  of 
50  KHz.  The  rest  of  the  nominal  input  specifications 
are  same  as  in  section  4.3.1,  except  that  for  one  case 
the  drift  tube  radius  has  been  halved  from  4  cm.  to  2  cm. 
The  dimensions  and  the  other  synchronous  parameters  are 

listed  in  table  2  and  table  3. 

It  is  seen  that  at  50  KHz  the  total  length  of  the 
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Table  2 

Structure  Dimensions  and  Synchronous  Parameters 


N 

Length  of  Nth. 
drift  tube 
(cm) 

Length  of 
Nth.  gap 
(cm) 

Transit  time 
factor 

Synchronous 

velocity 

(Km/Sec) 

Synchronous 

energy 

(Mv) 

1 

3.394 

.8485 

,2184 

4,243 

.3 

2 

3  c  572 

.8931 

.2441 

4,443 

,3284 

3 

3.743 

.9357 

22696 

4,657 

.3601 

4 

3,925 

.9811 

,2962 

4.883 

,3951 

5 

4,115 

1.0290 

.3135 

5,119 

.4336 

6 

4,312 

1.0780 

,3510 

5.366 

.4756 

7 

4,517 

1.1290 

.3785 

5.621 

.5212 

8 

4,727 

1.1820 

.4056 

5.883 

.5704 

9 

4,943 

1.2360 

.4322 

6.152 

.6231 

10 

5.162 

1.2910 

.4580 

6,426 

.6792 

11 

5.385 

1.3460 

,4829 

6.703 

.7387 

12 

5.610 

1.4020 

,5066 

6,984 

.8014 

13 

5,836 

1.4590 

,5293 

7.267 

.8673 

14 

6.063 

1.5160 

,5507 

7,551 

.9360 

15 

6,291 

1.5730 

.5710 

7.835 

1,0080 

16 

6.518 

1.6290 

.5900 

8.119 

1.0820 

17 

6,744 

1.6860 

,6080 

8.402 

1.1580 

18 

6.970 

1.7420 

,6248 

8,684 

1.2370 

19 

7.194 

1.7980 

.6406 

8.964 

1.3190 

20 

7.417 

1,8540 

.6554 

9.243 

1.4020 

21 

7.637 

1,9090 

,6692 

9.519 

1.4870 

22 

7.856 

1,9640 

,6822 

9,793 

1.5740 

23 

8.073 

2.0180 

,6944 

10.060 

1.6620 

24 

8,288 

2,0720 

,  7058 

10.330 

1.7530 

25 

8,501 

2.1250 

.7165 

10.600 

1.8440 

26 

8,712 

2,1780 

.  7265 

10.860 

1.9370 

27 

8,920 

2.2300 

.7360 

11,120 

2.0320 

28 

9.126 

2,2810 

.7448 

11.380 

2.1270 

29 

9.330 

2,3320 

,  7532 

11.640 

2.2240 

30 

9.531 

2.3830 

,7611 

11.890 

2.3220 

Total  length  of  the  structure  *  2.437  Meters . 

Note:  Operating  frequency  of  the  structure  =  50  KHz 

Other  input  specifications  are  same  as  in  section  4,3,1. 
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Table  3 

Structure  Dimensions  and  Synchronous  Parameters 


N  Length  of  Nth, 

Length  of 

Transit  time 

Synchronous 

Synchronc 

drift  tube 

Nth,  gap 

factor 

veloci ty 

energy 

(cm) 

(cm) 

(Km/Sec) 

(Mv) 

1 

3,394 

,8485 

,6360 

4,243 

,3 

2 

3,913 

,9784 

,6963 

4,827 

,3826 

3 

4,355 

1,0890 

,7388 

5,389 

,4731 

4 

4,781 

1,1950 

,7719 

5,923 

,5690 

5 

5,185 

1,2960 

,7978 

6,431 

,6693 

6 

5,570 

1,3930 

,8185 

6,915 

,7729 

7 

5,937 

1,4840 

,8353 

7,376 

o  8793 

8 

6,289 

1.5720 

,8493 

7,817 

,9878 

9 

6 , 626 

1,6570 

,8609 

8,241 

1,0980 

10 

6,951 

1,7380 

,8708 

8,648 

1,2100 

11 

7,263 

1,8160 

,8793 

9,040 

1,3230 

12 

7,566 

1,8910 

,8867 

9,419 

1,4370 

13 

7,858 

1,9650 

,8931 

9,786 

1,5520 

14 

8,142 

2,0350 

,8987 

10,140 

1,6680 

15 

8,417 

2,1040 

,9038 

10,490 

1,7850 

16 

8,685 

2,1710 

,9082 

10,820 

1,9030 

17 

8,946 

2.2370 

,9122 

11,150 

2,0210 

18 

9,200 

2,3000 

,9159 

11,470 

2,1390 

19 

9,449 

2,3620 

,9191 

11,780 

2,2580 

20 

9,691 

2,4230 

,9221 

12,080 

2.3770 

21 

9,929 

2,4820 

,9248 

12,380 

2,4970 

22 

10,160 

2,5400 

,9273 

12,670 

2,6170 

23 

10,390 

2,5970 

,9297 

12,960 

2,7380 

24 

10,610 

2,6530 

,9318 

13,240 

2,8590 

25 

10,830 

2,7080 

,9337 

13,510 

2.9800 

26 

11,050 

2,7620 

,9356 

13,780 

3,1010 

27 

11,260 

2,8140 

,9373 

14,050 

3.2230 

28 

11,470 

2,8660 

,9389 

14,310 

3.3440 

29 

11,670 

2,9170 

,9403 

14,560 

3.4660 

30 

11,870 

2,9680 

,9417 

14,810 

3,5880 

Total 

length  of  the 

structure 

-  3,137  Meters, 

Note : 

Inner  radius 

of  the  drift  tubes  -  2  cm 

Operating  frequency  of 

the  structure  - 

50  KHz 

Other  input  specifications  are  same  as 

in  section 

4,3,1 

- 
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structure  for  30  sections  is  almost  half  of  that  of  the 

structure  designed  for  30  KHz.  It  is  important  to  notice 

that  for  30  sections  the  energy  gain  for  the  structure 

designed  at  30  KHz  is  much  larger  than  that  for  the 

structure  designed  for  50  KHz,  in  the  case  where  the 

drift  tube  radius  has  been  held  constant.  The  reason 

for  this  behavious  is  that  the  ratio  a/LS  increases 

n 

and  hence  the  transit  time  factor  decreases,  as  may 
be  seen  from  equation  Al-10  of  Appendix  1.  On  the  other 
hand  it  may  also  be  seen  from  Table  3  that  if  the  drift 
tube  radius  for  the  50  KHz  structure  is  halved,  then 
the  transit  time  factor,  and  hence  the  overall  energy 
gain  is  greatly  increased  over  than  shown  in  Table  2. 

From  fig.,  4.8  it  is  seen  that  the  particle  injected 
on  the  axis  of  the  50  KHz  structure  at  a  given  slope 
collides  with  the  24th  drift  tube  wall  of  the  structure, 
if  there  are  no  quadrupoles .  If  the  particle  is  injected 
at  the  same  slope  into  the  structure  with  quadrupoles 
(voltage  between  quadrupoles  is  40  Kv)  it  strikes  the 
quadrupole  surface  of  the  21st  drift  tube.  The  reason 
is  that  although  the  transverse  excursions  of  the  particle 
are  less  in  the  structure  with  quadrupoles,  the  maximum 
permissible  transverse  excursion  is  decreased  because  of 
the  reduced  aperture.  Thus  the  net  gain  in  introducing 
the  quadrupoles  at  40  Kv  is  negative.  However,  the  gain 


- 
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EXCURSION  IN  VERTICAL  PLANE  EXCURSION  IN  HORIZONTAL  PLANE 

(In  Meter  X10)  (ln  Meter  X10) 

-0.30  -0.15  -0.00  0.15  0.30  _-0.30  -0.15  -0.00  0.15  0.30 
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Without  Quadruples 


With  Quadrupoles 


(Voltage  between  the  Quadrupoles=40  Kv) 


With  Quadrupoles 
(Voltage  between  the  Quadrupol es=100  Kv) 


13.00  19.00 

N 


25.00 


31.00 


Without  Quadrupoles 


With  Quadrupoles 


(Voltage  between  the  Quadrupoles=40  Kv) 


-  With  Quadrupoles 
(Voltage  between  the  Ouadrupoles=  100  Kv) 


- 1 - 1 — 

.00  7.00 


H - 1 - 1 - 1 - 1 - 1 - 1 

13.00  19.00  25.00  31.00 

N 


Note:  The  injection  conditions  of  the  particle  are  the 
same  as  for  the  trajectory  a  of  Fig.  4.2. 

FIG.  4.8  Transverse  Excursions  for 
Structure  With  and  Without 
Quadrupoles  Designed  at  50  KHz 
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in  introducing  quadrupoles  is  positive  if  the  quadrupole 
voltage  is  increased.  In  fig.  4.8  the  particle  traverses 
all  the  sections  without  interception  with  a  quadrupole 
voltage  of  100  Kv. 


4.3.7  A  Comparison  of  the  Equations  of  Motion  Derived 
in  This  Thesis  and  Those  used  by  Previous  Workers 


The  numerical  values  of  -Ar  ,  AZ  and  AE  /q  given 

n  n  n'  ^  ^ 

by  equation  2-8,  2-20  and  2-22  respectively  are  computed 
for  a  particle  injected  on  the  axis  of  the  accelerator 
at  a  certain  injection  slope.  The  structure  used  is 
the  one  designed  in  section  4.3.1.  The  above  equations 
are  new  in  the  sense  that  they  are  better  approximations 
to  the  particle  motion  in  an  accelerator  gap,  than  those 
equations  that  have  been  used  by  previous  workers  (ref¬ 
erences  6,  7,  8  and  14).  The  advantage  of  the  above 
equations  is  that  they  are  valid  for  any  radial  excursion 
and  that  the  transit  time  factor  is  correct  for  any 
nonsynchronous  velocity.  The  results  of  the  computation 

are  listed  in  table  4.  The  same  quantities  _A^nold 

AE 

/\z  and  11011:1  given  by  the  equations,  used  by  previous 

nold  q 

authors,  as  shown  below  are  computed  for  the  same  set 
of  injection  conditions  and  the  same  structure  as  above. 
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'  AWd  =  (q/m)  (n/LSn)  (Sin*n/vn)  (Vm)  (TSin(r)) - (4-i) 

AWd  =  (q/m>  <C0SW  (Vm)  (TSfn(r>) - (4-2) 

AEnold/cJ  =  (C0SV  (V  (Ts,n(r)) - (4-3) 

The  results  of  the  above  computation  are  given  in  table  5 

Upon  comparing  table  4  and  table  5  it  is  found 
that  the  results  of  the  computation  by  the  two  different 
sets  of  formulas  agree  very  well  for  the  first  few 
sections.  However,  in  later  sections,  when  the  particle 
no  longer  travels  near  the  axis  and  its  velocity  is  no 
longer  synchronous,  the  difference  between  the  two  sets 
of  formulas  becomes  important.  This  is  evident  from  the 
results  given  in  the  two  tables.  The  parameters  have 
been  computed  up  to  the  24th  section  because  the  particle 
is  intercepted  by  the  25th  drift  tube  wall. 
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Table  4 

Numerical  Values  of  the  Excursions,  and  Ar  .  AZ  and 
.  *  n'  n 

AEn/q  Given  by  the  Equations  2-8,  2-20  and  2-22 


N 

Arn 

(Meter/sec) 

^n 

(Km/sec) 

^Wv") 

X  (N) 

(cm) 

Y(N) 

(cm) 

1 

-.416 

.  4884 

.06907 

.0 

.0 

2 

-1.155 

.  4884 

.07702 

.01006 

.'01006 

3 

-1.823 

.  4784 

.08324 

.02111 

.  02111 

4 

-2.470 

.  4656 

.08842 

.03390 

.03390 

5 

-3.135 

.4513 

.09272 

.04923 

.04923 

6 

-3.844 

.  4367 

.09630 

.06785 

.06785 

7 

-4.619 

.4225 

.09930 

.09059 

.09059 

8 

-5.476 

.  4088 

.10180 

. 11830 

. 11830 

9 

-6.433 

.3959 

. 10400 

. 15200 

. 15200 

10 

-7.504 

.  3838 

.10590 

.19270 

. 19270 

11 

-8.705 

.3725 

. 10760 

.24150 

. 24150 

12 

-10.050 

.  3620 

. 10900 

. 30000 

.30000 

13 

-11.570 

.3522 

. 11030 

. 36950 

.36950 

14 

-13.260 

.3431 

. 11150 

.45170 

.45170 

15 

-15.170 

.  3347 

. 11260 

.54850 

.54850 

16 

-17.320 

.  3269 

. 11370 

.66210 

.66210 

17 

-19.740 

.3198 

.11470 

. 79480 

. 79480 

18 

-22.470 

.3133 

.11570 

.94930 

.94930 

19 

-25.570 

.  3073 

.11680 

1.12900 

1.12900 

20 

-29.090 

.  3020 

.11790 

1.33600 

1.33600 

21 

-33.130 

.  2974 

. 11920 

1.57600 

1.57600 

22 

-37.800 

.2934 

,12060 

1.85200 

1.85200 

23 

-43.240 

.2902 

. 12220 

2.17000 

2 . 17000 

24 

-49.660 

,2879 

. 12420 

2.53500 

2 . 53500 

Note:  The  particle  is  injected  at  an  angle  of  .01  degree  with 

the  axis  of  the  accelerator  designed  in  section  4.3.1. 
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Table  5 


Numerical  values  of  the  Excursions  and  Ar  ,  ,, 

nold 


AZ  ,  ,  and  AE  .  _/q 
nold  nold  ^ 

Given  by 

the  Equations 

• 

i — 1 

• 

■=3* 

2  and  4 . 3 

N 

• 

Ar  i  j 
nold 

AZ 

n 

AEnold/<3 

X  (N) 

Y(N) 

(Meter/sec) 

(Km/Sec) 

(Mv) 

(cm) 

(cih) 

1 

-.416 

.  4884 

.06907 

.0 

.0 

2 

-1.155 

.  4884 

.07702 

.01006 

.01006 

3 

-1.823 

.  4784 

.08324 

.02111 

.02111 

4 

-2.470 

.  4656 

.08842 

.03390 

.03390 

5 

-3.135 

.4513 

.09272 

.04923 

.04923 

6 

-3 . 845 

.  4367 

.09629 

.06785 

.06785 

7 

-4,620 

.  4225 

.09930 

.09059 

.09059 

8 

-5.478 

.40*88 

. 10180 

.11830 

. 11830 

9 

-6 . 436 

.  3959 

.10400 

. 15200 

. 15200 

10 

-7.509 

.3838 

.10590 

. 19270 

. 19270 

11 

-8.714 

.3725 

. 10760 

. 24160 

. 24160 

12 

-10.070 

.  3620 

. 10900 

. 30000 

.  30000 

13 

-11.590 

.3522 

. 11030 

. 36960 

. 36960 

14 

-13.300 

.  3431 

.11150 

. 45180 

.45180 

15 

-15.230 

.  3347 

. 11260 

.54880 

. 54880 

16 

-17.410 

.  3269 

.11360 

.66250 

.66250 

17 

-19 . 880 

.3198 

. 11460 

.79540 

. 79540 

18 

-22.690 

.3133 

.11560 

.95030 

.95030 

19 

-25.890 

.3074 

. 11670 

1.13000 

1.13000 

20 

-29 .580 

.  3021 

. 11780 

1.33900 

1.33900 

21 

-33.870 

.  2974 

. 11890 

1.58000 

1.58000 

22 

-38.910 

.2935 

. 12030 

1. 85800 

1.85800 

23 

-44.900 

.  2904 

. 12180 

2.17900 

2.17900 

24 

-52.1500 

.  2881 

. 12370 

2 . 54900 

2 . 54900 

Note : 

The  particle 

is  injected  at  an  angle 

of  .01 

degree  with 

the  axis  of  the  accelerator  designed  in  section  4.3.1. 
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Chapter  5 


CONCLUSION 


It  has  been  shown  that  a  Sloan-Lawrence  drift 
tube  structure  designed  for  a  particle  with  a  fixed 
nominal  charge  to  mass  ratio  accepts  particles  with  a 
wide  range  of  charge  to  mass  ratios,  both  higher  and 
lower  than  the  nominal  value.  The  range  of  acceptable 
charge  to  mass  ratio  for  a  particle  injected  at  certain 
slope  with  the  axis  of  the  accelerator  is  largest  if 

the  structure  contains  quadrupole  focusing  elements. 

for 

For  instance,/ a  structure  with  quadrupoles  and  designed 
for  a  charge  to  mass  ratio  of  30  Coulombs/Kg,  the 
acceptance  range  is  approximately  26.2  Coulombs/Kg  to 
83.5  Coulombs/Kg.  The  exact  acceptance  range  depends 
upon  the  particular  injection  conditions. 

The  energy  gain  of  a  particle  in  a  structure 
consisting  of  a  given  number  of  drift  tube  sections 
depends  critically  upon  the  field  geometry  via  the 
ratio  of  drift  tube  radius  to  section  length.  For 
instance,  for  a  30  drift  tube  section  structure  designed 
for  a  charge  to  mass  ratio  of  30  Coulombs/Kg  and  operat¬ 
ing  at  50  KHz,  the  energy  gain  of  a  particle  is  increased 
by  more  than  50%  if  the  drift  tube  radius  is  halved 
from  4  cm  to  2  cm. 
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The  inclusion  of  quadrupoles  increases  the 
acceptance  for  particles  which  are  injected  at  an 
angle  to  the  accelerator  axis .  For  a  structure 
designed  for  a  given  charge  to  mass  ratio  and  a  fixed 
quadrupole  voltage,  the  focusing  effect  of  the 
quadrupoles  is  best  at  lower  frequencies.  For  higher 
frequencies  the  section  lengths  and  consequently  the 
quadrupole  lengths,  are  shorter,  and  hence  the  focus¬ 
ing  effect  is  decreased. 

A  new  set  of  equations  to  describe  the  particle 
motion  in  the  accelerator  has  been  developed.  These 
equations  are  more  accurate  for  the  off  axis  and 
nonsynchronous  particles  than  those  equations  used  by 
previous  workers.  For  a  given  accelerator  structure 
and  a  given  set  of  injection  conditions,  the  two 
different  sets  of  equations  lead  to  results  which  differ 
by  as  much  as  5%. 
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APPENDIX  1 
TRANSIT  TIME  FACTOR 

In  the  Sloan-Lawrence  structure  the  "transit  time  factor"  for  a 
synchronous  particle  in  the  nth  section  is  defined  as  follows: 


LSn/2 


E 7(r,Z)  COS-™-  dZ 

L  L  On 


•LSn/2 


LSn/2 


Ez(r,Z)  dZ 


•LS  / 2 
r 


- - - - (Al-1) 


In  order  to  compute  T  (r)  ,  E  (r,Z)  is  evaluated  subject  to  the  folio- 

s  ^  n  z. 

wing  boundary  conditions: 


Ez(a,Z)  -  0 

=  vm/gn 


for  -LSn/2<z  <  -gn/2 
for  -gn/2<z<  gn/2 


gn/2<z  <  LSn/2 


(Al-2) 


Because  the  microparticle  accelerator  operates  at  a  very  low  fre¬ 
quency,  the  wave  length  'X  of  the  electromagnetic  field  is  very  much 
larger  than  the  length  of  a  gap»  Hence,  retardation  times  in  a  gap 
are  negligible  and  the  gap  field  satisfies  Laplace's  equation.  Thus 
the  scalar  potential  V(r,Z)  in  the  gap  is  given  by 
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2 

V  V(r9Z)~  0 

The  solution  of  the  above  equation  is 


V(r9Z)  =  A  I0(pr)+B  KQ(pr)j  C0S(pZ)+  D  SIN(pZ) 
where  A,  B,  C9  D  and  p  are  constants  and,  I  and  Kq 
Bessel  functions  of  first  kind. 


Since,  V(r,Z)  is  finite  at  r=0,  B  is  zero.  Thus  V(r 


are  the  modified 


,Z)  becomes 


V(r,Z)  -  |YcOS(pZ)  +  K^SIN(pZ)  I  (pr) 


where  K^nd  K*are  the  new  constants. 

It  is  now  assumed  that  the  structure  is  periodic  in  Z.  Although  the 
actual  period  is  approximately  LSn,  for  the  purpose  of  the  computation 
of  the  gap  fields  it  is  assumed  that  the  periodicity  is  (m)(LS  ),  where 
m  can  be  taken  to  be  any  even  integer.  The  later  assumption  has  little 
effect  on  the  actual  gap  field,  but  it  greatly  simplifies  the  comput¬ 
ations,  Thus, if  m  is  taken  as  2,  then  p  in  the  above  equation  can  be 
written  as 


KIT 


LSn 

where  K  =  1 ,  2 s  ------- 

The  foregoing  equation  now  becomes 

V(r,Z)  =  •{  k'cOS (— 3”)  +  K"siN(Sp-) 


LS. 


In(— )  - 
°  LSn 


(Al-3) 


'n  "  / 

V( a,Z)  as  specified  by  Al-2  is  now  expanded  in  a  Fourier  Series  as 


To  satisfy  the  above  boundary  condition  Al-3  is  written  as 


■ 
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00 


KIIZ 


V(r,Z)  =  g  J  Ak  COStj^)  +  B  SIN(^)U_V^LS 

n  n  i/tt 


Kllr, 


UtT") 


n  > 


y  V  J  f  Klla  \ 
LOn 


(Al-4) 


Since  the  electrical  centre  of  a  gap,  Z  =  0,  is  defined  as 


LS  /2 
n 


j  Ez(r,Z)  SIN(S|  )  dZ  =  0 


(Al-5) 


•ls  n 

n 


It  follows  that  V(r,Z)  is  an  odd  function  of  Z.  Thus  Al-4  becomes 


00 


Knr . 

I  77“ 

w/  „  0  LSn  ,Knz  . 

V(r,Z)  =  E  Bk  n  SIN(l$  ) 

K=1 

Hence, 


,  ,  3V(r,Z)  „  KH  wr  > 

Ez(r*Z)  =-  32 - ■  -E 


K=1 


I  (■^-) 

0  LSn  „„C,KHZ  , 

- COSfj-rr-  ) 

I  (KM)  LSn 

0  LS„ 


Substituting  Al-7  in  Al-1  one  obtains 


(Al-6) 


(Al-7) 
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•LSn/2 


I  (J®C) 

KE  °  LS" 

K=1  B|<  LSn  Jjm.) 


0  LS 
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COS(^)  COS ({2-  )  dZ 
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s  ,n 


LS/2 


-LSn/2 


Ez(r,Z)  dZ 
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Since  -  0  for  K  even,  one  obtains 
LSn/2 

-LV2 

Now,  assuming  that  the  particle  traverses  the  section  at  constant  r 

and  using  the  foregoing  result, T  (r)  becomes 

s  j  n 


B  COS  COS(S-)  dZ  =  0  for  K  f  1 

K  Lsn  LSn 


T 

s,n 


(r) 


-(Bi ) 


n 


1  LS. 


cos2  (5|  )  dZ 

Li)n 


n 


i  1 

VlsJ 


(B-i )  (g)  — JiT- 
l  ^  T  (M  ) 

0\s  ’ 

LOn 


V 


(Al-8) 


m 


B  is  found  out  from  the  boundary  condition  Al-2  as  follows: 
From  Al-2  and  Al-6  one  obtains 


LS  /2 
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-LSn/2 


,nz 

LS 


gn/2 


n 


—  C0S(5|  )  dZ 
9n  LSn 


00  i/TT  kiizx 

-Z  B  (§■  )  COSfj^-)  dZ 
K=1  K  Lin  n 


-LSn/2 


Upon  integration  the  last  two  integrals  become 


■  ni  o  r. 
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Vm  4  LSn 

(—)  ( 

9n 


,ngn 

n  ^  SIN^2LSn^  =  "  Bin 


-  bx  =  ~  (^)  sin(5?J1  ) 

1  n2  gn  v2LSn; 


(Al-9) 


From  Al-8  and  Al-9  the  transit  time  factor  can  now  be  written  as 


T 

s  ,n 


(r) 


s,,(lr„>  'o<5„> 


Ei  Vg ' 

2LS  n 


(Al-10) 
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APPENDIX  2 

INTEGRATION  OF  EQUATIONS  2-7  OF  PAGE  9  AND  2-19  OF 

PAGE  14  TO  COMPUTE  A r  AND  AZ 

n  n 

To  perform  the  integration  of  equation  2-7  one  proceeds  as  follows: 
The  radial  component  of  the  electric  field  is  computed  from  Al-6 
of  appendix  1  as 


E  (r,Z)  = 
r 


9V(r,Z) 

3r 


.KJlr 

VrsT1 


„  vu  ,  )  - - -  SIN(j^p-) 

K=1  K  LSn  j  (MS.  )  Lsn 


=  -  E  (B  )  (tt  ) 


(A2-1) 


‘0VLS 


n 


KHrv  . 


where  I  is  the  modified  Bessel  function  of  first  order. 

1  LSn 

Thus,  the  first  integral  on  the  R.H.S.  of  equation  2-7  can  be 
written  as 


LS  JZ 


nz 


E  (r,Z)  SIN(jj  )  dZ  = 


■LS  / 2 
n 


KTIr 

LSn/2  MlS^ 

o°  Kn 

1=1  K  LSn  ,  (M) 
-LSn/2  ’0lLSn) 


SIN(]^)SIN(2|  )  dZ 
Li>n 


Assuming  that  the  particle  traverses  the  gap  at  constant  r  the 
above  equation  becomes 
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LV2 


-LS  /2 
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Er(r,Z)  SIN(I!|  )  dZ  = 


LS 


I  (—  ) 
^1  ISfi 


I  (M  ) 
0  LS 


Substituting  the  value  of  from  Al-9  and  using  the  equation  Al-10 
the  foregoing  equation  becomes 


LS  /2 
n 


f  Er(r,Z)  SIN(nZ/LSn)  dZ  =  Vffl  Ts>n(0)  I^JIr/LS^  . (A2-2) 


-LS  /2 
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This  can  be  written  as 


LSn/2 


-LSn/2 


Er(r,Z)  SIN(nZ/LSn)  dZ  =  R  V 


m 


(A2-3) 


where  R  =  T$ ^(0)  I^nr/LS^ - - - — (A2-4) 

To  find  the  second  integral  of  equation  2-7  one  proceeds  as 
follows : 

Differentiating  both  sides  of  equation  A2-3  w.r.t.  LSn  one  obtains 


where  Er(r,LSn/2)  and  Ep( r LSn/2)  are  the  values  of  Er  at  the  two 
limits  Z  =  ±  LSn/2.  Since  a  very  little  error  is  introduced  in 
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assuming  that  the  limits  of  integration  iLS^/2  can  be  extended  infinite¬ 
ly  on  both  sides,  one  can  assume  in  the  foregoing  equation  to  be 
zero.  Thus  the  above  equation  becomes 
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d_ 

dLSn 


|Er(r,Z)  SIN(JIZ/LSn)|  dZ  = 
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dR 

dLS, 


-LSn/2 


Upon  differentiation  one  obtains 
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f  (HZ/LSn)  Er(r,z)  C0S(nZ/LSn)  dZ  = 


■  V  LS  — 
vm  L^n  dLS 


(A2-5) 
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-LSn/2 


dR 

dLS. 


of  A2-4  is  given  by 


dR  _  £ 
dL$n"LSn 


siN(52a) 

ZLbn 


I  ) 
A0  LSn 


fn9n  v 

2LSn 


(“) 

LS 

L  n 


Ilr  Ttr 

1 1 (ltr/LSn)  +  rj-  I2(^-) 
_  Lbn  c 


/  TEr  »  T  /Hi  \ 
‘  a  ll  LS^  !l'LSn 


nr 


>+ 


It  (—  ) 
1  LS. 


I  (—  ) 
0  LS„ 


n  /9n,2 


.(Hi)"  C0S(^^')“—  SIN  (Hill  ) 
LSn  2  2LS'  2  2LS„ 


rng, 


(5!a ) 

2LS„ 


(A2-6) 


110 


Tq  find  the  integral  of  equation  2-19  one  can  write  the  equation  Al-1 
of  appendix  1  as 


LSn/2 


J 


E,(r,Z)  C0S(nZ/LS  )  dZ  =  Vm  T  (r) 
l  n  m  s  ,n 


( A2-7 ) 


-LSn/2 


where  T  (r)  is  given  by  Al-10  of  appendix  1. 

b  j  FI 

Now,  differentiating  both  sides  of  A2-6  w.r.t.  LSn  one  obtains 
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where  *  is  given  by 
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APPENDIX  3 
PHASE  OSCILLATION 


The  synchronous  phase  changes  by  2n  in  going  from  one  section 
to  the  next.  The  change  in  phase  for  a  nonsynchronous  particle  moving 
from  the  n*'*1  to  the  (n+l)^  section  with  velocity  vn  is  given  by 


n  +  cot  =  n  +- 

n 


2n  LSn(l+kn) 


ToV 


s  ,n 


=  n  (2+kn)  - - (A3-1) 


th 


The  nonsynchronous  phase  at  the  centre  of  the  n  section  is 


♦n  =  V  ^(2+kP) 


(A3-2) 


where  $  is  the  phase  of  the  electric  field  at  the  instant  the  non- 
o 

synchronous  particle  is  injected. 

The  synchronous  phase  at  the  centre  of  the  n*'*1  section  is  given  by 


cj)  =  6  +  2nH 
s,n  ys 


where  <j>  is  the  phase  of  the  electric  field  at  the  instant  the  syn- 
s 

chronous  particle  is  injected.  Thus,  the  oscillation  of  nonsynchronous 
phase  about  the  synchronous  phase  is 
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=  (*  -  *  )+  n 

0  S 


n 

Z  At 

p=l 


+  kp)  -  2n  k 
) 


< 


— (A3-4) 


